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ABSTRACT 

The  perturbation  method  proposed  by  Professor  Hui  is 
described.   The  method  gives  exact  solutions  for  the  perturbed 
flow  over  both  sides  of  a  flat  plate  which  is  oscillating 
with  small  amplitude  and  frequency  at  large  angles  of  attack 
in  steady  supersonic/hypersonic  inviscid  flow  provided  that 
the  shock  remains  attached.   Using  the  strip  theory  concepts 
these  solutions  are  extended  to  study  the  dynamic  stability 
in  pitch  of  a  flat,  periodically  oscillating  wing  or  arbi- 
trary planform  shape,  at  large  angles  of  attack.   Finally, 
Hui ' s  perturbation  method  is  extended  to  include  the  effects 
of  upstream  disturbances  on  a  stationary  wedge. 
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I.   INTRODUCTION 

This  thesis  deals  with  oscillating  bodies  in  inviscid, 
steady  or  unsteady,  supersonic/hypersonic  flow.   Its  sub- 
ject matter  may  be  divided  into  three  main  topics. 

The  first  topic  is  covered  in  Sections  II  and  III. A. 1,2 
and  constitutes  the  background  material  for  the  other  two. 
In  Section  II  the  Eulerian  governing  equations  and  boundary 
conditions  for  unsteady  flow  are  formulated  and  the  two- 
dimensional  shock  and  expansion  steady  flow  results  are  given. 
The  basic  elements  of  the  linearized  potential  flow  theory 
are  also  included  in  this  section.   In  Sections  III .A. 1,2  a 
perturbation  method  proposed  by  Professor  Hui  is  presented 
[Ref.  6,7,8].    The  method  uses  as  a  basis  the  assumption 
that  the  unsteady  flow  over  an  oscillating  flat  plate,  with 
attached  shock  waves  at  an  arbitrary  angle  of  attack,  is  a 
small  perturbation  from  the  steady  reference  flow  and,  for 
small  amplitudes  of  periodic  oscillations,  it  gives  closed 
form  solutions  for  the  flowfield  quantities  in  the  disturbed 
flow  regions. 

The  second  topic  is  covered  in  Sections  III.A.3,B,C 
[Ref.  5].   In  Section  III. A. 3  the  closed  form  solutions  over 
the  upper  and  lower  sides  of  the  oscillating  plate  are  combined 


Professor  of  Applied  Mathematics,  University  of  Waterloo, 
Ontario,  Canada.   NPS,  Department  of  Aeronautics,  Visiting 
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and  the  strip  theory  approximation  is  employed  to  calculate 
the  necessary  quantities  in  the  disturbed  flow  regions  over 
a  flat  wing  of  arbitrary  planform  shape.   The  in-pitch  stability 
derivatives  for  the  three-dimensional  wing  are  then  obtained 
in  closed  form  and  the  results  are  compared  with  other  exist- 
ing theories.   In  Section  III.C  a  comparison  of  these  results 
with  linearized  potential  flow  theory  results  is  included 
while,  in  Section  III.B,  the  fundamentals  of  the  linearized 
theory,  as  applied  to  three-dimensional  wings,  are  presented. 
Finally,  the  third  topic  is  covered  in  Section  IV.   Based 
on  the  same  perturbation  method  this  topic  introduces  upstream 
unsteadiness  in  the  flow  and  its  effects  on  stationary  bodies 
are  obtained.   The  upstream  unsteadiness  is  of  a  fairly  general 
periodic  form  to  give  the  solution  in  the  flowfield  generated 
by  a  formation  of  bodies  provided  that  the  body  originating 
shocks  are  not  crossing  and  the  expansion  fan  regions  are 
not  overlapping.   An  extension  to  the  case  of  oscillating 
bodies  is  readily  possible. 
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II.   UNSTEADY  INVISCID  FLOW  THEORY 

To  describe  the  fluid  motion  two  methods  are  available: 
the  Lagrangian  method  and  the  Eulerian  method.   In  both 
methods  the  fluid  is  regarded  as  a  continuum,  i.e.,  its 
matter  is  assumed  to  be  continuously  distributed. 

In  the  first  method  the  fluid  is  assumed  to  be  divided 
into  infinitesimally  small  regions  called  fluid  elements  or 
fluid  particles.   The  so-called  particle  point  of  view  is 
then  adopted  and  description  of  the  fate  of  each  individual 
fluid  particle  is  sought .   To  determine  the  unknowns  asso- 
ciated with  each  fluid  particle,  e.g.,  its  position  coor- 
dinates, density  etc.,  a  system  of  equations  is  set  up  by 
applying  to  each  fluid  particle  natural  laws  such  as  Newton's 
second  law  of  motion  and  conservation  of  mass  and  energy. 
These  equations  are  known  as  the  Lagrangian  equations  of 
fluid  motion. 

Although  the  Lagrangian  description  appears  to  be  a 
natural  way  to  approach  the  problem  of  fluid  motion,  the 
Eulerian  description  is  preferred  in  general  since  it  gives 
more  insight  into  the  problem,  it  is  much  simpler  and  in 
most  cases  one  is  not  interested  in  the  fate  of  each  individual 
fluid  element  but  rather  in  the  properties  of  the  fluid  at  a 
certain  point  of  the  flow  field  at  a  certain  time. 

In  the  Eulerian  method  attention  is  focused  on  the  vari- 
ous points  of  the  space  filled  by  the  flowing  fluid  and  a 
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description  of  what  is  happening  at  each  of  these  points, 
in  terms  of  quantities  of  interest  such  as  pressure,  density 
and  velocity  is  sought.   The  flow  quantities  of  interest 
are  assumed  to  be  functions  of  space  and  time,  i.e.,  to  be 
scalar  or  vector  fields.   Thus  in  the  so-called  field  point 
of  view  adopted  here  the  fluid  flow  is  characterized  by  the 
fields  of  velocity,  pressure,  density  and  so  on  and  a  fluid 
element  occupying  a  certain  point  at  a  certain  time  assumes 
for  its  properties  the  values  that  are  appropriate  to  that 
point  at  that  time.   To  solve  for  these  fields  a  system  of 
equations  is  again  set  up  by  using,  as  before,  natural  laws 
such  as  Newton's  second  law  of  motion  and  conservation  of 
mass  and  energy.   These  equations  are  known  as  the  Eulerian 
equations  of  fluid  motion. 

Throughout  this  thesis  the  Eulerian  approach  is  used. 

In  the  next  subsection  the  major  steps  in  deriving  the 
Eulerian  equations  of  fluid  motion  are  indicated  and  the 
equations  are  presented  in  the  form  in  which  they,  will  be 
used  later  on. 

A.   UNSTEADY  EULER  EQUATIONS 

The  Eulerian  equations  of  fluid  motion,  called  hereafter 
simply  equations  of  motion,  may  be  set  up  either  in  differen- 
tial form  or  in  integral  form.   Furthermore  they  may  be 
developed  either  from  the  point  of  view  of  a  certain  fluid 
region  that  contains  the  same  fluid  elements  for  all  times 
(control  mass  approach)  or  from  the  point  of  view  of  a  fixed 
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volume  in  space  through  which  different  fluid  elements  flow 
through  (control  volume  approach) . 

In  what  follows,  the  equations  of  motion  are  set  up  in 
the  differential  form  from  the  point  of  view  of  an  infinitesi- 
mal fluid  region. 

In  this  derivation  we  will  then  naturally  be  involved  with 
the  calculation  of  the  time  rate  of  change  of  flow  quantities 
as  we  follow  the  fluid  element  around,  the  so-called  material 
derivatives  of  quantities. 

A  physical  interpretation  of  the  material  derivative-  and 
its  components  is  briefly  included. 

1.   Material  Derivative 

Consider  a  fixed  coordinate  system  and  a  fluid  ele- 
ment situated  at  point  r  at  some  time  t.    Let  Q(r,t)  denote 
some  fluid  property  Q  of  interest  (density,  velocity,  etc.) 
associated  with  the  point  r  at  time  t.   The  fluid  element 
situated  there  (see  Fig.  1.1)  will  assume  for  its  corresponding 
property  Q  the  value  Q(r,t).   In  a  short  time  interval  At 
the  element  moves  through  a  distance  As  =  VAt  where  V  is 
its  velocity  at  r  and  t.   The  element  will  then  assume  for 
Q  the  value  appropriate  to  its  new  position  r  +VAt  at  time 
t+At.   If  we  denote  this  value  as  Q(r+VAt,  t  +  At)  ,  the 
change  of  Q  in  the  time  interval  At  is 

AQ   =  Q(r  +VAt,t  +  At)  -Q(r,t) 


Barred  quantities  denote  vector  quantities 
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and  the  rate  of  change  of  Q  following  the  element  around, 
usually  denoted  by  p^   is 

DQ      ,.    Q(r+VAt,t+At)  -  Q(r,t) 
Dt   "   At*0  At 

Expanding  Q (r+VAt, t+At)  in  a  Taylor  series  we  get 


Q(r+VAt,t+At)   =   Q(r,t)  +  (|§)    At  +  C^-4)    At2  +  ... 

dZ   r,t      8t  r,t 

2 

+  (|§L   VAt  +  c^-4)     (vAt)2  +  ... 

s  r,t       8s   r,t 

where  s  denotes  distance  in  the  direction  of  the  velocity 
V  at  point  r  and  time  t. 

Using  the  Taylor  series  expansion  the  rate  of  change  of 

n   DQ  , 

Q,  jr^-  becomes 


where  the  derivatives  are  evaluated  at  r  and  t. 

The  total  rate  of  change  of  any  property  Q  is  thus  com- 
posed of  two  parts . 

To  see  the  physical  interpretation  of  each  of  these  terms 
consider  a  flow  field  which  is  at  any  instant  spacewise  uni- 
form but  varies  from  instant  to  instant  and  a  second  flow 
field  which  is  steady  but  not  uniform  spacewise.   Consider 
also  a  fluid  element  in  these  fields  which  in  a  small  time 
interval  5t  is  moving  from  position  r  in  the  flow  field  to 
position  r+VSt. 
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As=VAt 


poth 


(r*Ar)  Qt    (UAt) 


Fig.  1.1.   Illustration  of  local,  convective  and 
material  derivatives. 


The  change  of  a  property  Q  of  the  fluid  element  moving 

30 
in  the  first  field  described  above  is  -r^-6t  and  the  rate  at 

d  t 

which  the  property  Q  is  changing  locally  at  the  point  r  is 
•£-£■.   This  local  rate  of  change  is  known  as  the  local  deriva- 
tive and  is  the  first  component  of  the  total  rate  of  change 
in  equation  (1.1). 

The  change  of  a  property  Q  of  the  fluid  element  moving 

30 
in  the  second  flow  field  described  above  is  (^-=-)V5t.   This 

dS 

change  which  is  called  convective  change  is  necessary  since 
the  element  has  to  have  an  appropriate  value  for  its  property 
Q  at  its  new  position  r+V5t.   The  rate  of  change  of  the 
property  is  (■*-=•)  V  and  is  known  as  the  convective  derivative. 

d  S 

It  is  the  second  component  of  the  total  rate  of  change  in 
equation  (1.1) . 


17 


The  sum  of  the  local  derivative  and  the  convective 
derivative  as  given  by  equation  (1.1)  is  known  as  the  total 
or  substantial  or  material  derivative.   The  last  term  is 
probably  more  descriptive  since  the  derivative  is  constructed 
following  a  certain  material  element  around.   This  term  will 
be  used  hereafter. 

Recalling  that  -!—   represents  the  derivative  of  Q  with 

a  S 

respect  to  distance  in  the  direction  of  the  velocity  V  we  can 
write  equation  (1.1)  in  the  following  form 


8|  -  H  ♦  <ev.gradQ)V 


where  e   is  the  unit  vector  in  the  direction  of  V  or  simply 
§|  =  fl  +  V-gradQ  tl.2) 


where  Q  can  be  a  scalar  or  a  vector  quantity. 

If  Q  is  vector  quantity,  i.e.,  Q  =  A  the  convective 
derivative  V«gradA  can  be  expanded  using  the  formula 

V-gradA  =  ^-[grad(V-A)  -  V  x  curlA  -  A  x  curl  V 

-  curl(VxA)  +V(divA)  -A-(divV)]         (1.3) 

2 .   Equations  of  Motion 

We  will  now  set  up  the  basic  equations  that  govern 
the  unsteady  motion  of  an  inviscid,  compressible  fluid.   We 
initially  regard  as  unknowns  the  velocity  field  V(r, t) ,  the 
pressure  field  p(r,t)  and  the  density  field  p(r,t).   We 
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want  to  establish  relationships  between  these  fields  by 
applying  to  a  certain  fluid  element  the  basic  laws  of  nature; 
Newton's  second  law  of  motion,  law  of  conservation  of  mass, 
law  of  conservation  of  energy, 
a.   Momentum  Equation 

Let  us  consider  an  inf initesimally  small  fluid 
element  situated  at  position  r  at  time  t  (Fig.  1.2) .   If 
V  and  p  are  the  velocity  and  density  of  the  element  at  r  and 
t  and  if  6v  denotes  the  volume  of  the  element,  then  the  mass 
and  the  momentum  of  the  element  are  p6v  and  p6vV  respectively. 
Let  us  also  denote  by  F  the  total  force  acting  on  the  element 
at  time  t. 

Newton's  second  law  of  motion  which  is  applicable 
to  any  mechanical  system  states  that  "at  any  instant,  the 
rate  of  change  of  momentum  of  a  system  is  equal  to  the  resultant 
of  all  forces  that  are  acting  on  the  system  at  that  instant". 
By  applying  the  above  law  to  the  fluid  element  considered 
and  by  noting  that  the  rate  of  change  of  momentum  of  the  ele- 
ment is  simply  the  material  derivative  of  the  momentum  we 
get 


§£(p6W)   =   F 


Since  the  mass  of  the  element  p6v  remains  constant  the 
above  equation  becomes 


P6v  §|  =  F  (1.3.1) 
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-pndS 


Figure  1.2.   Fluid  Element 

The  total  force  F  is  the  resultant  of  the  so- 
called  surface  forces  and  body  forces . 

The  body  forces  are  forces  that  act  throughout  the 
body  of  the  fluid  such  as  the  gravity  force.   These  forces 
will  be  assumed  to  be  small  and  will  be  neglected  in  this 
thesis . 

The  surface  forces  on  the  other  hand  are  internal 
forces  in  the  nature  of  actions  and  reactions  across  the  sur- 
fact  that  separates  the  fluid  element  from  its  neighboring 
fluid  elements.   For  a  frictionless  or  inviscid  fluid  the 
surface  forces  are  simply  pressure  forces  that  act  normal  to 
the  surface  of  the  fluid  element.   Their  resultant  is 

-  <£$   pndS 
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where  p  is  the  pressure  of  the  element,  5S  its  total  surface 
and  n  the  outward  unit  normal  to  the  surface. 

According  to  the  integral  definition  of  the  gradi' 
ent  of  a  scalar  function  the  resultant  of  the  pressure  forces 
which  is  the  total  force  F  can  be  written 


F   =   -Cft>pndS   =   -6v  gradp 


Equation  (1.3.1)  then  becomes 


pz;  =   -gradp  (1.4) 


This  equation  of  motion  is  one  of  the  fundamental  equations 
of  fluid  dynamics  and  is  called  Euler's  Equation.   It  repre- 
sents a  system  of  three  scalar  equations  for  the  five  un- 
knowns— the  pressure,  the  density  and  the  three  scalar  com- 
ponents of  velocity. 

b.   Mass  Equation 

Consider  a  fluid  element  situated  at  point  r  at 
time  t  with  volume  5v  (Fig.  1.3) .   The  mass  of  the  element 
is  p6v.   The  law  of  conservation  of  mass  states  that  "the 
mass  of  any  fluid  element  remains  constant  as  it  moves  about" 
even  though,  in  general,  its  shape,  volume  and  density  may 
change . 

Applying  the  law  of  conservation  of  mass  to  the 
element  considered  is  thus  equivalent  to  setting  the  rate 
of  change  of  mass  equal  to  zero  or  setting  the  material 
derivative  of  mass  equal  to  zero.   This  gives 
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or 


g^(p<5v)   =   0 


X    DP         D    X 

6V  Dt  +  P  Dt  6V   = 


(1.5) 


The  material  derivative  of  the  volume  of  the  element  may  be 
expressed  in  terms  of  the  velocity  field.   Let  the  surface 
of  the  element  at  time  t  be  6S  and  let  it  in  a  time  inter- 
val 6t  grow  and  become  6S,  as  shown  in  (Fig.  1.3) .   The 
change  in  volume  of  the  element  is  equal  to  the  volume  swept 
by  the  surface  of  the  element  during  the  time  St.   If  n  is 
the  outward  normal  to  the  original  surface  the  net  volume  swept 
outward  by  6S  in  time  5t  is  given  by 


Figure  1.3.   Change  in  volume  of  a  fluid  element 
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jfc^VStnds 

where  V  is  the  velocity. 

The  material  derivative  of  the  volume  of  the 
element  is  then  if  we  also  employ  the  definition  of  the 
divergence  of  a  vector 


jjr-  6v   =   (£1 


Equation  (1.5)  then  becomes 


|£  +  p  div  V   =   0  (1.6) 


This  equation  is  known  as  the  equation  of  conservation  of 
mass,  or  simply,  the  equation  of  mass,  or  the  equation  of 
continuity.   It  is  a  relation  between  the  velocity  and  den- 
sity fields  only.   Since  it  does  not  involve  any  dynamical 
quantities  (such  as  pressures  or  forces)  it  is  a  kinematical 
relation. 

c.   Energy  Equation 

The  law  of  conservation  of  energy  expresses  the 
balance  of  energy  exchanges  that  take  place  between  a  system 
and  its  surroundings.   A  fluid  in  motion  may  be  regarded 
as  a  thermodynamic  system  characterized  by  the  usual  thermo- 
dynamic variables  such  as  entropy,  internal  energy,  etc. 
We  will  assume  that  the  fluid  is  non-heat  conducting  and 
also  that  for  a  fluid  element  the  only  possible  energy  ex- 
change process  is  work  done  by  the  surface  forces  and  body 
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forces.   The  law  of  conservation  of  energy  applied  to  a 
fluid  element  may  be  expressed  as  follows. 

The  rate  of  increase  of  energy  E  of  a  fluid 
element  =  the  rate  of  work  W,  done  by  the  surface  forces  and 
the  rate  of  work  W2  done  by  the  body  forces.   Symbolically 

5§  =  wx  +  w2  (1.7) 

We  will  assume,  as  before,  that  the  body  forces  can  be 
neglected  and  that  the  fluid  is  inviscid.   The  only  possible 
energy  exchange  is  then  work  done  by  the  pressure  forces. 
The  rate  of  this  work  is 

W,   =   -  t/pndS-V  =   -  <rtpV-ndS 

1        6S  5S 

where  6S  as  before,  is  the  surface  of  the  element,  n  is  the 
outward  unit  normal  and  p,  V  are  the  pressure  and  velocity 
of  the  fluid  element. 

According  to  the  definition  of  the  divergence  of 
a  vector  we  can  write 

W,   =   -  (C^pV-ndS   =   -6v  div(pV) 

6S 

On  the  other  hand  the  energy  E  of  the  fluid  element 
is  the  sum  of  its  kinetic  energy  and  internal  energy.   We 
specify  the  internal  energy  of  the  fluid  by  the  scalar  field 
e(r,t)  which  denotes  the  internal  energy  per  unit  mass  at 
point  r  and  time  t.   Then  since  the  kinetic  energy  per  unit 
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v2 

mass  is  -j-   the  total  energy  of  the  fluid  element  with  mass 
p6v  is 

E   =   p6v(e  +  ^-) 


Equation  (1.7)  can  then  be  expressed  as 

2 

P  §^(e  +  ~)   =  -  div(pV)  CI. 8) 

Equation  (1.8)  is  referred  to  as  the  equation  of  conserva- 
tion of  energy  or  simply  the  energy  equation.   An  alternative 
form  of  this  equation  is 


p  ||  =   -  p  div  V  (1.9) 


which  can  be  found  by  subtracting  from  (1.8)  the  so-called 
equation  of  mechanical  energy 

2 
P  Dt(~2")   =   "  v*9radP 

The  equation  of  mechanical  energy  is  formed  by  multiplying 
both  sides  of  Euler's  equation  (1.4)  by  V. 

The  energy  equation  (1.8)  has  introduced  the 
internal  energy  of  the  fluid  as  an  additional  unknown  in  the 
formulation  of  the  governing  equations.   The  list  of  unknowns 
includes  the  three  scalar  components  of  velocity,  and  also 
the  pressure,  density  and  internal  energy  of  the  fluid,  while 
there  are  five  equations  available. 
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At  this  point  we  assume  we  are  dealing  with  a 
perfect  gas  and  introduce  the  equation  of  state  for  a  per- 
fect gas 

p   =   P  R  T  (1.10) 

as  a  sixth  relation  between  the  unknowns. 

We  specify  the  temperature  T  of  the  fluid  as  a 
scalar  field  T(r,t)  and  since  we  are  dealing  with  a  perfect 
gas  express  the  internal  energy  e  of  the  fluid  by  the  relation 


e   =   Cv  T  (1.11) 


where  C   is  the  specific  heat  at  constant  volume  of  the  gas. 

From  equations  (1.9),  (1.10)  and  (1.11)  we  can 
find  that  the  energy  equation  for  a  perfect  gas  may  be 
written  in  the  following  form  which  will  be  used  hereafter 

}fc<i)   =   0  (1.12) 

Summarizing  we  state  that  the  basic  equations  that  govern 
the  unsteady  motion  of  a  non-heat  conducting,  inviscid,  per- 
fect gas  with  constant  specific  heats  are  equations  (1.4), 
(1.6)  and  (1.12).   These  equations  are  rewritten  below  for 
easy  future  reference. 


(continuity)       f£  +  V-(pV)   =   0  (la) 

(momentum)         j^  +   V-gradV  +  ^2-  =  0            (lb) 

(energy)           ft^  +  7*V(^}   =  °             Clcl 
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3.   Boundary  Conditions 

Physical  conditions  that  should  be  satisfied  on 
given  boundaries  of  the  fluid  are  known  as  boundary  condi- 
tions.  There  are  several  types  of  boundaries  and  consequently 
there  are  various  possibilities  for  the  boundary  conditions . 

We  will  consider  two  types  of  boundaries  which  are 
of  main  importance  in  this  thesis,  CD  "the  solid-fluid 
boundary"  where  the  fluid  is  bounded  by  a  solid  surface  and 
(2)  "the  fluid-shock-fluid  boundary"  where  two  regions  of  the 
same  fluid  in  different  states  of  motion  are  separated  by  a 
flow  discontinuity.   The  possibility  of  an  infinitely  weak 
discontinuity  will  not  be  excluded. 

The  nature  and  number  of  the  boundary  conditions  depend 
also  on  the  form  of  the  differential  equations  that  govern 
the  motion  of  the  fluid.   In  this  sense  there  are  differences 
between  the  boundary  conditions  for  a  viscous  fluid  or  an 
inviscid  fluid.   In  the  following  the  conditions  for  an  invis- 
cid  fluid  are  considered. 

a.   Conditions  at  a  Solid-fluid  Boundary 

We  assume  that  the  fluid  is  bounded  by  an  imper- 
meable solid  wall  and  require  that  no  fluid  should  cross  the 
solid  surface.  Since  the  surface  itself  may  be  in  motion  we 
denote  by  V  the  velocity  of  the  fluid  and  by  V  the  velocity 
of  the  surface.  The  relative  velocity  between  the  fluid  and 
the  surface  is  V-V  .  Let  the  equation  of  the  surface  be 
given  by  Str,t)  =0.  A  unit  normal  to  the  surface  is  then 
given  by 
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—      gradS 
tjgradS 


and  the  component  of  the  relative  velocity  normal  to  the 
surface  is  given  by 


tV-Vs).n  =   tlgAdFltV-gradS-V^gradS] 


Now  assume  that  an  observer  moves  with  the  surface  particles 
that  compose  the  solid  surface.   The  observer  cannot  observe 
any  change  in  the  function  S(r,t)  considered  as  a  scalar 
field.   This  means  that  the  total  rate  of  change  of  S(r,t), 
following  a  particle  of  the  surface  around,  is  zero,  i.e., 


||+Vs-gradS      =      0 


The  component  of  the  relative  velocity  normal  to  the  sur- 
face then  becomes 


(V-Vn  '  ±|gradF|(ff  +  V-gradS)  (1.13) 


The  condition  of  impermeability  of  the  solid  surface  is 


(V-V  ) -n   =   0 
s 


or 


g|  -  ||  +  V  gradS   =   0    At  S(r,t)  =  0    (1.14) 


If  the  solid-fluid  boundary  is  formed  by  the  surface  of  a 
stationary  rigid  solid  the  above  equation  reduces  to 
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V-gradS   =   0    At   S  (r)   =   0  (1.15) 

It  is  pointed  out  here  that  the  condition  for- 
mulated above  states  that  at  each  point  of  the  solid-fluid 
boundary  the  normal  to  the  surface  component  of  the  relative 
velocity  between  the  fluid  and  the  solid  must  vanish.   Thus 
for  an  inviscid  fluid  nothing  can  be  said  about  the  tangen- 
tial to  the  surface  component  of  the  relative  velocity  which 
may  or  may  not  be  zero.   In  short  the  so-called  no-slip  con- 
dition does  not  apply  to  an  inviscid  fluid. 

b.   Conditions  at  the  Fluid-Shock-Fluid  Boundary 

Consider  one  dimensional  adiabatic  constant-area 
flow  of  a  perfect  gas  through  a  discontinuity  (Figure  1.4a)  . 
Assume  that  the  flow  quantities  in  regions  1  and  2  are  con- 
stant throughout  the  regions.   The  equations  of  continuity, 
momentum  and  energy  between  cross  sections   1   and   2   give 
[Ref.  1:  pp.  55,56] 


P1U1   =   P2U2 


Pl  +  P1U1   =   P2  +  P2U2 


2  2 

ui     .,   Pi      u0     ^   p 


2    y-1  P-l      2    y-1  p2 

These  equations  hold  as  long  as  sections   1   and   2   are 
chosen  outside  the  discontinuity  region.   The  discontinuity 
region  may  be  assumed  to  be  vanishingly  thin  and  sections   1 
and   2  may  be  brought  arbitrarily  close  together.   In  this 
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Figure  1.4a.   Change  of  conditions  across  a  normal 
shock  in  a  constant  area  duct. 
Steady  flow. 


Region 


S(r,t) 


Figure  1.4b. 


Change  of  conditions  across  an  arbitrary 
ciscontinuity  surface  S(r,t).   Unsteady  flow 
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case  the  requirement  of  a  constant  area  duct  is  dropped  and 
the  equations  apply  locally  across  the  discontinuity.   Con- 
sidering the  discontinuity  as  a  boundary  separating  regions 
1  and  2  we  may  refer  to  equations  above  as  boundary  conditions 
at  the  discontinuity  and  write  symbolically 

[pu]   =   0  CI. 16) 

At  x  =  x 

2  ° 

[pu   +  p]   =   0  (1.17) 

tiT  +  ^l?]  -  °  (1'18) 

where  the  square  brackets  denote  the  change  in  the  enclosed 
quantity  across  the  discontinuity. 

A  generalization  of  the  simple  one-dimensional 
flow  problem  considered  leads  to  the  well  known  Rankine- 
Hugoniot  conditions  in  the  form  that  will  be  used  in  this 
thesis.   Consider  two  regions  in  space  separated  by  a  surface 
S(r,t).   Assume  that  adiabatic  flow  of  an  inviscid,  perfect 
gas  is  established  from  region  1  to  region  2  (Figure  1.4b) . 
Assume  also  that  the  surface  S(r,t)  represents  a  discontinuity 
In  general  the  equation  of  the  discontinuity  is  not  known 
a  priori  but  will  be  found  as  part  of  the  solution  of  the 
flow  problem.   Thus  the  boundary  is  a  so-called  free  boundary. 

Let  the  flow  quantities  in  region  1  be  V, (rft), 
p^Fjt),  p1(r,t)  and  in  region  2  V2(r,t),  p2(r,t),  p2(r,t). 
The  boundary  conditions  (1.16),  (1.17)  and  (1.18)  should  apply 
locally  at  any  point  of  the  boundary  S(r,t)  provided  that  the 


31 


velocities  u,,  u   are  replaced  by  the  normal  components  of 

velocities  V   ,  V     For  a  boundary  moving  with  velocity 

nl   n2 
V  the  normal  to  the  boundary  components  of  the  relative 

velocities  between  the  fluid  and  the  free  boundary  are  given 

by   formula    (1.13) 


V 
n. 


±lgradSl(?f  +  VgradS) 


V  =      — ■ — L__(||.  +  v    -gradS) 

n2  ± |gradS |    3t  2   * 

From  equations  (1.16),  (1.17)  and  (1.18)  using  the  same 
notation  we  get: 

At  S  =  0 
(continuity)   [p  (||  +  V«VS) ]   =   °  (1.19) 

(normal        [p (||  +  V- VS) 2  +  p (VS) 2]   =   0  (1.20) 

momentum) 

(energy)     [|(||  +  V- VS)  2  +  -Ij.  £( VS)  2]   =   0  (1.21) 


where  the  square  brackets  denote,  as  before,  the  change  in 
the  enclosed  quantity  across  the  discontinuity  and  the  symbol 
"V"  stands  for  "gradient". 

The  conservation  of  tangential  momentum  is  not 
expressed  by  any  of  the  equations  above.   To  find  the  tan- 
gential momentum  equation  we  require  that  the  velocity  com- 
ponent tangent  to'  the  discontinuity  be  continuous.   The 
tangential  velocity  components  are  given  by 
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V    =   (n  x  v  )  xn 


V    =   (n  x  v  )  xn 
^2  ^ 


where  n  is  the  outward  unit  normal  to  the  surface 

n  =  VS  =  {S  ,S  ,S  }  in  Cartesian  coordinates.   The  condition 
x  y   z 

of  conservation  of  tangential  momentum  then  becomes  if  the 
square  bracket  notation  is  used: 

(tangential    [u]   _   [v]   _   [w]  Q 

momentum)     S       S       S 

x       y       z 

where  [u] ,  [v] ,  [w]  denote  the  change  of  x,  y  and  z  components 
of  velocity  across  the  discontinuity. 

The  equation  above  imposes  two  scalar  conditions 
at  the  discontinuity,  as  one  should  expect  from  physical 
considerations.   An  alternative  form  of  this  equation  is  the 
following: 


S1     ^>VS,   -   °     At   S   =   0  (1.22) 


Equations  (1.19),  (1.2)),  (1.21)  and  (1.22)  constitute  the 
complete  set  of  the  Rankine-Hugoniot  conditions.   They  are 
symmetrical  and  therefore  remain  unchanged  if  the  brackets 
are  taken  to  denote  the  change  upstream  rather  than  downstream 
through  the  discontinuity.   A  definite  sense  of  flow  direction 
is  provided  by  the  second  law  of  thermodynamics ,  which  requires 
that  the  entropy  shall  not  decrease  across  a  discontinuity. 
The  change  of  entropy  is  given  by  [Ref.  1:  p.  60]: 
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R  Pj_  P]_ 

and  the  requirement  stated  gives  the  following  condition 
which  should  accompany  the  Rankine-Hugoniot  conditions 

(^™Vf        t-rl   >   0    at  S   =   0  (1.23) 

thermo . )  Y    — 

P 

It  should  be  noted  that  so  far  we  have  avoided  using  the 
term  shock  instead  of  discontinuity,  though  shocks  are  the 
only  possible  physical  discontinuities.   This  was  done  on 
purpose  since  we  intend  to  use  the  Rankine-Hugoniot  condi- 
tions across  hypothetical  expansion  fronts  (negative  or  expan- 
sion shocks)  through  an  iterative  procedure  so  that  in  the 
limit  condition  (1.23)  should  not  be  violated. 
4 .   Shock-Expansion  Flows 

Unsteady  flows  with  shock  waves  or  expansion  waves 
are  considered  in  this  thesis .   These  unsteady  flows  will  be 
solved  by  first  finding  the  corresponding  steady  flow  solution 
and  then  using  it  as  a  reference  flow  in  calculating  the  un- 
steady perturbation  flow. 

Since  a  number  of  exact  steady  flow  conditions  are 
already  available,  they  will  be  utilized  as  reference  flows 
in  finding  the  governing  equations  and  boundary  conditions 
of  the  corresponding  unsteady  flows.   They  include  the  super- 
sonic uniform  wedge  flow  and  the  Prandtl -Meyer  expansion  flow. 
The  results  for  these  steady  flows  are  stated  below  for  easy 
reference. 
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a.   Wedge  Flow 

Consider  steady  uniform  supersonic  flow  past  a 

symmetrical  two-dimensional  wedge  with  semi-vertex  angle  9 

(Figure  1.5).   The  wedge  is  assumed  stationary.   Oblique 

shocks  will  be  formed  at  angles  S  measured  from  the  free 

stream  direction.   By  conservation  of  momentum  the  tangential 

component  of  velocity  is  continuous  across  the  shock  so  that 

V\   =  V\  .   Then  V   and  V   are  related  by  the  normal  shock 
t]_    t2         n1  n2 

relations.   Since  V   =  V,  sin  8  and  V   =  V_  sine})  the  normal 

nx    1  n2    2 

shock  relations  can  be  used  directly  with  M.  replaced  by 

M  sin  3  and  M  replaced  by  MLsin<|>.   The  resulting  relations 


for  the  oblique  shocks  are 


2     2 

2yM.  sin  S  -  (y-1) 


Y+l 


M,.P,.p, 


V, 


Figure  1.5. 


Steady  supersonic  flow  past  a  stationary 
wedge 
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2  2  (Y-l)M2sin23  +2 

M_    sin    <J>      = 


2yM2sin23  -  (y-D 


P2  (y+l)M2sin23 


pl  2  +  (y-l)M2sin23 


p2  tan   B 


p,  tan  <$> 


From  the  last  two  relations  eliminating  —  and  recalling 
that  <J)  =  6-9  we  get 


p2 


pl 


2    2 

M'sinl  -1 

tan  9   =   cot  3  — n ^ 5 5 (1.24) 

Ili  M,  -  (MTsin^B  -1) 


One  way  to  solve  this  equation,  i.e.,  to  find  3  for  given 

M,  and  0  is  to  express  it  as  a  cubic  equation  in  x  =  cot  3  and 

select  the  appropriate  positive  root  corresponding  physically 

to  the  weak   shock  wave  [Ref.  2:  pp.  452-453].   The  following 

equivalent  equation  was  used  in  this  thesis  for  numerical 

calculations 


(1  +^M2)  x3  -  (M2-l)cot9  x2  +  (1  +^  M2)x +cot  9   =   0 


where  x  =  tan  3 . 

For  attached  shock  waves  this  equation  gives  three 
real  roots  for  3  and  the  middle  one  is  the  one  corresponding 
to  the  weak  shock  wave. 
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A  maximum  deflection  angle   for  given  M.  (beyond 
which  the  shock  becomes  detached  )  can  be  found  from  equation 
(1.24)  by  differentiating  it  with  respect  to  3  and  equating 
to  zero. 

b.   Prandtl-Meyer  Flow 

Consider  steady,  two  dimensional,  uniform  super- 
sonic flow  over  a  convex  corner  (Figure  1.6a).   A  turn  of 
the  flow  through  a  single  oblique  expansion  wave  is  not 
possible  since  this  would  lead  to  a  decrease  in  entropy. 
The  flow  expands  isen tropically  through  an  infinite  number 
of  centered  straight  Mach  lines  that  form  the  so-called 
Prandtl-Meyer  expansion  fan.   Thus  upstream  of  the  ray  OB 
where  9=9   the  flow  is  uniform  with  Mach  number  M   and 

00 

downstream  of  the  ray  OC  where  9=9,  the  flow  is  also  uni- 
form with  Mach  number  M, .   For  angles  9  such  that  9   <  9  <  9, 

1  oo    _         _       i 

the   flow   field  has    the   same  properties    along   any   ray 
9   =   constant.      The  polar   coordinate    system  shown  has   been 
chosen  so   that   r   =  a   +  5-  +  P  (M    )    where  P  (M)    is    the   Prandtl- 
Meyer   function 


P(M)      =      j  tan   1[X(M2-1)1/2]    -    tan"1 (M2-l) 1/2 

y-1  1/2 
with  A  =  (-h^r)     and  y    the  ratio  of  specific  heats  of  the 
Y+l 

gas .   The  Mach  number  M,  is  given  by 


For  the  wedge  at  an  angle  of  attack  a  the  flow  deflec- 
tion angles  are  9+a  and  9-ct  for  the  lower  and  upper  sides 
respectively. 
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M„.p»  -p. 


u. 


t/JJ/f///  /  /  /  ?'JJ/tlill 


Figure  1.6a.   Steady  supersonic  flow  over  a  convex  corner 


m^rrrm 


actual    surface 


hypothetical   surface 


stagnant    fluid  region 


Figure  1.6b.  Expansion  over  a  corner  with  a  >  I (M) 


max 
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P  (M,  )   =   P  (M  )  +  a 

1         °° 


Along  any  ray  9  in  the  expansion  fan  the  Mach  number  pressure 
and  density  are  given  by  the  following  relations. 


-  j-   tan~1[A(M2-l)1/2] 


E(M  )   ,.       , . 
P   =   p  r !L.]Y/(Y-1) 


E(M  )  ,  ,.       ,. 

p   =  p  r^^_]V(Y-l) 


00 


E(M) 


where 


y-1   2 
E(M)   =   1  +  ■!-—   M 


The  r-  and  6-velocity  components  in  the  expansion  fan  are 
given  by 


where : 


V   =   c  sin 

r 


V\   =   Ac  cos 


=  0  y1+  2    i 


Y-l  M2 

1        J^OO 


It  can  be  seen  that  the  Prandtl-Meyer  function 
P  (M)  defined  above  has  a  maximum  value  for  M  •*■   °° .   This  value 
is 
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P(M)max   =   J(/(Y+13/ty-l)-l) 


For  convex  corner  angles  a  >  P(M)max  -  POyM  the  streamlines  behind 
the  expansion  fan  behave  as  though  the  flow  occurred  over  an 

expansion  of  P (M)     and  for  an  inviscid  fluid  a  region  of 

c  max  ' 

stagnant  fluid  lies  between  the  hypothetical  position  of  the 
body  as  sensed  by  the  flow  and  the  actual  position  of  the 
body  (Figure  1.6b).   Whenever  this  occurred  in  numerical 
calculations  performed,  the  influence  of  the  flow  on  the  sur- 
face of  the  corner  was  assumed  negligible. 

B.   LINEARIZED  POTENTIAL  EQUATION 

The  so-called  linearized  theory  of  supersonic  flow  builds 
up  the  flow  produced  by  the  motion  of  a  body  by  superposi- 
tion of  small  disturbances  such  as  those  produced  by  a  moving 
sound  source.   One  can  develop  in  this  way  relatively  simple 
methods  for  the  computation  of  velocity  and  pressure  distri- 
butions in  the  field. 

In  the  case  of  vortex- free  flow,  the  equations  of  motion 
can  be  reduced  to  equations  analogous  to  the  wave  equation. 
The  coordinate  parallel  to  the  direction  of  the  main  flow 
plays  the  role  of  the  time  coordinate.   Hence  the  methods 
of  finding  solutions  of  the  wave  equation  can  be  used. 


We  mean  that  in  the  equivalent  case  of  a  flat  plate 
at  an  angle  of  attack  a  the  pressure  on  its  upper  surface 
was  assumed  zero. 
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The  linearized  theory  however  has  serious  limitations. 
First,  it  gives  only  a  first  approximation  since  all  devia- 
tions from  the  uniform  parallel  flow  are  considered  small 
and  therefore  additive.   This  is  justified  for  very  thin  or 
slender  bodies  at  small  angles  of  attack  only.   Second,  there 
are  speed  ranges  in  which  the  linearization  of  the  equation 
of  motion  even  for  small  disturbances  is  not  justified.   For 
the  linearized  theory  to  be  valid  the  following  two  conditions 
must  be  met. 

(a)  The  perturbation  velocities  must  be  small  in  compari- 
son to  both  the  main  stream  velocity  and  the  velocity  of  sound 
This  condition  excludes  the  case  of  very  high  velocities 
since  if  the  mean  stream  velocity  is  several  times  larger 
than  the  sound  velocity,  disturbances  which  are  small  rela- 
tive to  the  mean  stream  velocity  may  be  of  the  same  order  of 
magnitude  as  the  sound  velocity.   This  speed  range  is  called 
the  hypersonic  range . 

(b)  The  perturbation  velocities  must  be  small  in  comparison 
to  the  difference  of  the  main  stream  velocity  and  the  sound 
velocity.   This  condition  excludes  the  range  near  M  =  1,  the 
so-called  transonic  range. 

In  spite  of  the  limitations  described  above  the  lineari- 
zation of  the  equations  of  inviscid,  compressible  fluids 
proved  to  be  of  excellent  use  in  developing  approximate  solu- 
tions in  the  supersonic  range. 

There  are  three  general  methods  used  in  the  linearized 
theory  of  supersonic  aerodynamics. 
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(a)  The  method  of  fundamental  solutions  or  sources.   This 
method  is  based  on  the  superposition  of  fundamental  solutions 
of  the  linearized  hyperbolic  equation  for  small  perturba- 
tions of  a  uniform  supersonic  flow.   In  formulating  wing 
problems  this  method  uses  sources  and  doublets  located  in 
the  plane  of  the  wing  and  the  strength  of  these  singularities 
is  determined  so  that  the  boundary  conditions  applicable  to 
the  wing  planform  and  shape  are  satisfied. 

(b)  The  methods  of  acoustic  analogy  and  operational  cal- 
culus.  In  these  methods  the  solution  of  the  hyperbolic  equa- 
tion is  expressed  by  means  of  Fourier  and  Laplace  integrals 
respectively.   The  second  method  is  better  adapted  to  super- 
sonic flow  problems  since  the  Laplace  integrals  exclude  the 
possibility  of  upstream  travelling  signals  while  in  the  case 
of  Fourier  integrals  one  has  to  impose  additional  conditions 
to  secure  that  this  possibility  is  excluded. 

(c)  The  method  of  conical  flows.  This  method  is  based  on 
conical  flows,  i.e.,  flows  for  which  the  velocity  components 
at  points  lying  on  a  straight  line  drawn  from  a  point  chosen 
as  vertex  are  independent  of  the  distance  from  the  vertex. 

In  this  method  the  solution  of  a  hyperbolic  equation  in  three 
variables  is  reduced  to  the  solution  of  Laplace's  equation  or 
wave  equation  in  two  variables  and  the  existing  methods  of 
conformal  transformations  and  the  theory  of  functions  of  com- 
plex variables  can  be  employed. 

Methods  of  higher  approximations,  i.e.,  methods  which 
lead  from  the  simple  case  of  the  linearized  solution  toward 

o 
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the  exact  one  in  successive  steps,  thus  extending  the  range 
of  satisfactory  approximation,  are  required  whenever  the 
perturbation  velocities  are  not  small  compared  to  the  main 
stream  velocity. 

In  this  section  the  concept  of  irrotationality  will  be 
introduced  and  then  the  basic  assumptions  and  steps  followed 
in  deriving  the  linearized  potential  equation  and  the  appli- 
cable boundary  conditions  for  the  general  flow  problem  of  a 
body  in  supersonic  flow  will  be  described.   In  the  next  sec- 
tion the  linearized  wing  problem  will  be  described. 

1.   Irrotational  Flows 

By  potential  flow  we  mean  that  the  velocity  V  is 
derivable  from  a  scalar  velocity  potential  <j> ,  i.e.,  v  =  gradf 
On  the  other  hand  the  vorticity  or  rotation  uT  in  a  fluid  is 
defined  as  oT  =  curl  V  and  the  flow  is  called  irrotational 
if  a)  =  0  or  equivalently  if 


_9y_8u      3wav_2u_lw_     Q 
3x    3y      dy         3z      3z  "  9x 


where  u,  v,  w  are  the  x-,  y-  and  z-  components  of  V. 

Physically  the  irrotationality  of  the  fluid  means 
that  the  fluid  particles  have  zero  moment  of  momentum  about 
their  own  center-of— gravity' axes  or  simply  that  they  remain 
parallel  to  themselves  as  they  move  around. 

The  condition  of  irrotationality  of  the  flow  is  a 
necessary  and  sufficient  one  for  the  assumption  of  a  poten- 
tial flow  since  the  mathematical  identities  curl  grade)?  =  0 
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and  div  curl  V  =  0  show  that  the  velocity  V  can  then  be 
put  in  the  form  V  =  grad<t> . 

The  irrotationality  throughout  a  flow  field  can  be 
proved  by  applying  the  theorems  of  Kelvin  and  Stokes. 
Stokes1  theorem  states  that  the  area  sum  of  the  rotation  over 
a  given  area  is  equal  to  the  integral  of  the  velocity  around 
a  curve  bounding  the  area.   Formally, 

/  iJ  dA  =   /  curl  V-dA  =  <p  V-dT 
A 

The  line  integral  is  called  circulation  and  is  denoted  by 
T.   Thus  r  =  0  V-dl  =  JuT-dA. 

As  a  consequence  of  Stokes*  theorem,  co  =  0  if  the 
circulation  r  vanishes  for  all  paths  wholly  within  a  simply 
connected  flow  region. 

Kelvin's  theorem  on  the  other  hand  states  that  the 

circulation  r  about  any  contour  always  composed  of  the  same 

fluid  particles  (i.e.,  a  fluid  line)  is  constant  in  an  inviscid 

fluid  with  only  conservative  or  irrotational  body  forces. 

For  an  inviscid  fluid  it  states  that  =rr-  =  -(fi-£  and  it  reduces 

Dt     '    p 

DP 
to  =rr  =  0  when  there  is  a  simple  relation  connecting  p  and 

p.   Physically  it  means  that  circulation  r  about  any  line 

contour  remains  constant  in  time  as  we  move  along  with  the 

fluid. 

It  is  a  consequence  of  both  theorems  that  initially 

irrotational  flows  originating  in  a  reservoir  under  uniform 

stagnation  conditions  or  from  straight  parallel  streamlines 
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will  remain  irrotational  throughout  the  flow  field  at  all 
times  if  there  are  no  shock  waves. 

2.   Linearized  Potential  Flow  Equation 

In  view  of  the  above  reasoning  we  assume  that  the  whole 
flow  field  is  irrotational  and  set  V  =  grad$.   As  a  result 
the  three  unknown  velocity  components  u,  v,  w  are  expressed 
in  terms  of  the  scalar  field  $  by 


u   =   $  ,     v=<f>,     w  =   (l> 
x  y  z 


By  using  relation  (1.3)  the  momentum  equation  (1.4) 
may  be  written  as 


_         2 

p  (^r-  +  grad  -* —  V  *  curl  V)   =   -gradp 


or  since  curl  V  =  0  and  v  =  grad$  we  may  write 


,<3<I>   V2    f  dp, 
grad(_  +  __  +  /  _£)   = 


By   integration  we   get 


%+4+J±      -      F(t. 


If  we  define  $  =  <f>-/F(t)dt  this  equation  becomes 


By  differentiating  this  relation  with  respect  to  time 
and  by  taking  its  gradient  we  get  the  following  two  relations/ 
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if  the  additional  relation  a   =  -*£•  for  the  velocity  of 


sound  is  used. 


3tA3t    2  '      p  3t 
-  grad(—  +  -j)      =     —   gradp 


We  now  introduce  continuity  equation  (1.6)  which,  since 

_    2 

div  V  =  V  $,  may  be  written  as 


1  |£L  +  1  gradp  +  v2$  =  0 

p  at   p  * 


Introducing  in  this  equation  the  last  two  relations  from 
momentum  equation  we  get 

2       2  2 

V2$  -  A.[L-|  +  |X_  +  v  grad(^-)]   =   0  (1.25) 

a  dt 

Equation  (1.25)  is  the  exact  non-linear  differential  equation 
to  be  satisfied  by  the  velocity  potential  $  for  an  unsteady, 
inviscid,  irrotational  flow.   Because  of  its  strong  non- 
linearity,  solutions  have  been  found  in  very  few  special 
cases.   Thus  the  small  disturbance  concept  is  introduced  which 
leads  to  linearization  of  the  equation. 

We  assume  that  the  velocity  vector  V  differs  only 
slightly  in  direction  and  magnitude  from  the  free  stream 
velocity  U  ,  taken  along  the  x-axis,  and  we  define  a  distur- 
bance velocity  potential  ¥  obtained  from  the  total  velocity 
potential  $  by  separating  out  the  contribution  of  the  uniform 
flow.   In  this  way  we  have 
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$   =   Ux  +  ? 

The  local  velocity  components  are  then  given  by 

u  =   U   +  u',     v  =   v'  ,    w  =  w1 

where  u' ,  v' ,  w'  are  the  perturbation  or  disturbance  veloci- 
ties which  can  be  found  from 


u'   =   fx,     V   =   f  , 


We  regard  all  disturbance  velocities  small  in  compari- 
son with  U^,  a  and  U^-a  and  all  pressure  and  density  changes 
small  in  comparison  with  main  stream  pressure  and  density. 
We  also  assume  that  the  small  quantities  change  gradually 
in  all  directions  and  that  time  variations  are  not  too  rapid. 

We  now  return  to  equation  (1.25)  and  assume  that 
the  linear  terms  in  Laplace's  operator  are  of  the  same  order 
of  magnitude.   Substituting  the  velocity  V  in  this  equation 
and  retaining  first  order  terms  only  we  get 

2  2  2 

V2Y  -  -i[^-2r  +  2U  Wr   +  U2  i-1]   =   0  (1.26) 

where  the  second  term  in  parenthesis  is  found  from 


a  a  a 

and  the  third  term  in  parenthesis  is  found  from 
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-  4r[$2$     ]     =    -  -i-[(u  +uf)2y     ]     =    -  JL[u2^     ] 

2   x  xx         2    °°      xx         2   °°  xx 

a  a  a 

To  complete  the  linearization  of  (1.26)  we  should  dispose 

of  the  factor  —j. 

a  2 

For  steady  flow  a  may  be  eliminated  by  using  the 

following  relation  which  is  a  consequence  of  the  first  law 

of  thermodynamics 


2      Y_1  T72  .     . 

a   +  J-^—  V   =   constant 


For  unsteady  flow  let  a  slender  body  move  in  the  flow 
field  and  denote  the  velocity  of  sound  far  upstream  in  the 

undisturbed  field  by  the  constant  a^.   In  the  vicinity  of  the 

2  .  . 

body  a  will  then  be  a  variable  depending  on  position  and 

can  be  represented  by  a  sum  of  terms  of  the  form 

2     2       2 
a  -    a   +  (Aa)   +  ... 

00 


As  a  first  approximation  we  set  a  =  a   and  introduce  this 
value  in  equation  (1.26)  getting 

2  2  2 

V2y    -   -L [L4  +  2U      |-I-+U2    i-X]       =      0  (1.27) 

2    -.2  0°    8x9t        oo-2 

a      8t  3x 

00 

Uoo 

An  alternative  form  of  this  equation  with  M  =  —  is 

"*  oo        a 

oo 

?  2M  1 

(1  -M   )¥       +¥       +  V 4f    .    -  -=r  Y,..,.      =      0  (1.28) 

oo'    xx        yy         zz        a^       xt      ^2      tt 

00 

Equation  (1.28)  is  the  linearized  unsteady  potential  equa- 
tion and  is  used  as  the  basic  equation  in  most  aerodynamic 
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analyses.   Because  of  the  assumptions  made  and  terms  retained 
this  equation  is  valid  for  unsteady,  inviscid,  irrotational 
flows  that  are  purely  subsonic  or  purely  supersonic  and  is 
limited  to  small  disturbances  only. 

3 .   Linearized  Boundary  Conditions 

To  completely  specify  the  mathematical  problem  that 
describes  the  flow,  the  following  boundary  conditions  need 
generally  be  prescribed  [Ref.  3:  pp.  1.27-1.29]. 

a)  Surface  Boundary  Condition:   The  wing  surface  is 
impenetrable  to  the  medium. 

b)  Edge  Conditions:  Enough  viscosity  remains  in  the 
inviscid  fluid  to  determine  the  flow  pattern  near  sharp 
edges . 

c)  Wake  Conditions:   The  free  vorticity  shed  from  the 
trailing  edge  must  have  a  circulation  which  vanishes  together 
with  the  bound  circulation.   It  is  furthermore  assumed  that 
the  shed  wake  is  a  continuous  sheet  of  discontinuity  which 

is  coplanar  with  the  wing  projection  in  the  direction  of 
flight.   Edge  effects  and  rolling  up  of  the  sheet  are 
disregarded. 

d)  Conditions  at  Infinity:   A  state  of  uniform  flow  must 
be  prescribed  at  inifinity.   In  addition  the  Sommerfeld 
radiation  condition  requires  waves  to  propagate  away  from 
sources  of  disturbance  toward  infinity. 

e)  Other  Conditions:   As  the  most  important  additional 
condition,  the  requirement  that  proper  account  be  taken  of 
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zones  of  influence  and  action  at  supersonic  flow  velocities, 
is  mentioned. 

The  first  of  the  above  conditions,  namely,  the  con- 
dition of  impermeability,  as  applied  to  a  wing  or  airfoil, 
is  considered  next.   The  procedure  can  be  easily  extended  to 
cover  other  bodies  of  interest  such  as  slender  fuselages,  etc 

Consider  a  wing  fixed  relative  to  a  Cartesian  co- 
ordinate system  so  that  it  lies  close  to  the  xy-plane  (Figure 
1.7).   Assume  that  the  wing  is  submerged  in  an  infinite  mass 
of  fluid  moving  with  velocity  U^  in  the  positive  x-axis. 
Let  the  upper  and  lower  surfaces  of  the  wing  be  expressed  by 
equations 


S      z  -  z  (x,v,t)   =   0 
u  u   - 


S£   =   z  -  z^(x,y,t)   =   0 


The  condition  of  impermeability  of  these  surfaces, 
by  recalling  equation  (1.14),  requires  that 

9z      3z      8z 
w   =   TF  +  U  "3F  +  v  ~W     for   z  =  V   (x'y)  in  Ra 

dzl  dzl  dzl 

w  "  TF  +  u  TT  +  v  "ST  for  z  -  V   (x'^  in  Ra 

where  u,  v,  w  are  the  components  of  velocity  V  and  R   is 
the  portion  of  the  xy-plane  covered  by  the  projection  of  the 
planform. 

These  are  exact,  non-linear  equations.   To  linearize 
them  assume,  as  before,  that  the  disturbance  velocities  are 
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u. 


z  =  z 


mathematical    plane 
surface 


Figure  1.7.   Linearized  boundary  conditions  for  a  wing 


small  compared  to  the  free  stream  velocity  U   and  also  that 

3z    3z 

the  slopes  -*— -,  -w~,  etc.,  are  very  small  compared  to  unity 
o  x    dy 

Then  retaining  first  order  terms  only  we  get 


3z 


w  = 


w   = 


u 


3z 


3t 

if* 

3t 


+  U 


+  U 


u 


3x 

if* 

3x 


for    z  =  z  ,   (x,y)   in   R 
u      J        a 


for    z  =  z    ,       (x,y)   in   Ra 


Since  z   and  z.    are  small  compared  to  the  wing  chord 
we  may,  as  a  further  step,  replace  the  actual  wing  with  an 
infinitesimally  thick  surface  of  discontinuities  in  u,  v,  w 
and  pressure  p.   With  this  mathematical  plane  surface  located 
on  the  xy-plane,  we  may  expand  w  in  Maclaurin  series  about 
its  values  just  above  and  below  the  xy-plane 
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w(x,y,zu,t)   =  w(x,y,0+,t)  + z^   9w(x^Y'°  &1    + 


w(x,yfZjl,t)   =  w(x,y,0~,t)  +z   8w(x^f°  /t}  + 


Using  the  same  arguments  as  before  the  higher  order 
product  terms  can  be  neglected  and  the  impermeability  or 
flow  tangency  conditions  take  the  following  linear  forms 

Sz        3z 
w(x,y,0  ,t)   =   -—*■  +  U^  -j£  (x,y)   in   R&      (1.29a) 

9Z£       9Z£ 
w(x,y,0  ,t)   =  ~  +   U   -—■     (xfy)   in   R      (1.29b) 

a  X.  °°   ox  a 
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III.   OSCILLATING  WINGS  OF  GENERAL  PLANFORM 
IN  SUPERSONIC/HYPERSONIC  FLOW 


Consider  a  uniform  (spacewise  and  timewise)  supersonic 
or  hypersonic  flow  of  an  inviscid  perfect  gas  with  constant 
specific  heats  past  a  flat  wing  of  a  general  planform  shape 
at  an  angle  of  attack  a.      Assume  that  the  wing  is  performing 
a  small  amplitude  slow  pitching  oscillation. 

The  problem  considered  is  to  find  the  unsteady  flow 
quantities  in  disturbed  regions  over  the  wing  and  thus  its 
stiffness  and  damping  derivatives . 

The  governing  equations  of  motion  are  given  by  equations 
(la-lc) ,  restated  below. 

||  +  V-(pV)   =   0  (2.1a) 

U  +  V-VV  +  ¥£■  =   o  (2.1b) 

dt  P 

!_  (-£.)  +  V-V(-£-)   =   0  (2.1c) 

at  y       y 

p'       p' 

where  p,  p,  V  and  y    are  the  pressure,  density,  velocity 
and  adiabatic  exponent  of  the  gas. 

The  flow  tangency  condition  to  be  satisfied  at  the  sur- 
face of  the  body  is  given  by  equation  (1.14) ,  restated  below 


If-  +  V-VS   =   0     At   S   =   0  (2.2) 

d  t 


where  S(r,t)  =  0  is  the  equation  of  the  body  surface 
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The  fluid-shock-fluid  boundary  conditions  to  be  satis- 
fied across  the  shock  are  given  by  equations  (1.19-1.22), 
restated  below. 

[p(||  +  V-VG)]   =   °       At  G   =   0  (2.3a) 

[p(||  +  V-VG)2  +  p(VG)2]   =  0  (2.3b) 


[|(||  +  V-VG)2  +  ^  £(VG)2]   =   0  (2.3c) 

[VxVG]   =   0  (2.3d) 

where  G(r,t)  =  0  is  the  equation  of  the  unknown  shock  shape 
and  the  square  brackets  denote  the  change  in  the  enclosed 
quantities  across  the  shock. 

Equations  (2.1-2.3)  are  nonlinear.   The  nonlinearity  of 
the  governing  equations  and  boundary  conditions  along  with 
the  existence  of  a  shock  with  an  unknown  shape,  contribute 
to  the  complexity  of  the  problem  considered. 

For  low  supersonic  Mach  numbers  and  very  low  angles  of 
attack  shock  waves  can  be  replaced  by  Mach  waves  and  the 
linearized  supersonic  potential  flow  theory  can  be  employed. 
The  problem  can  then  be  solved,  at  least  for  certain  groups 
of  planform  shapes,  by  fairly  general  methods.   The  fundamen- 
tals of  the  linearized  supersonic  potential  flow  theory  as 
applied  to  a  three-dimensional  oscillating  wing  are  presented 
in  Section  III.B. 

For  high  angles  of  attack  and/or  Mach  numbers  the  shock 
waves  become  strong  and  the  linearized  theory  cannot  be  used. 
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To  overcome  the  difficulties  encountered  in  this  case,  namely 
the  existence  of  the  shock  and  the  nonlinearity  of  the  equa- 
tions, Professor  Hui  proposed  the  use  of  a  perturbation  method 
in  which  the  unsteady  flow  field  is  regarded  as  a  small 
perturbation  to  some  reference  steady  flow.   Thus  the  unsteady 
flow  problem  is  to  be  solved  by  first  finding  the  corresponding 
steady  flow  solution  and  then  using  it  as  a  reference  flow 
in  calculating  the  unsteady  perturbation  flow.   The  solution 
of  the  three-dimensional  wing  by  this  method  will  be  presented 
in  Section  III. A. 

Finally,  in  Section  III.C  results  for  the  stability  deriva- 
tives are  presented.   A  comparison  with  linearized  potential 
flow  theory  results  is  also  included. 

A.   PROFESSOR  HUI'S  THEORY 

In  this  section  the  problem  of  dynamic  stability  of  a 
flat  wing  of  a  general  planform  shape  at  arbitrary  angles  of 
attack  in  steady  supersonic/hypersonic  flow  is  considered. 
The  wing  is  assumed  to  be  oscillating  in  pitch  with  small 
amplitude  and  frequency  and  the  bow  shock  be  attached  to  the 
body  at  all  times . 

The  problem  is  covered  in  [Ref.  5]  and  only  the  basic 
steps  will  be  included  here,  in  Section  III. A. 3.   Its  solu- 
tion is  based  on  the  assumption  of  an  inviscid  perfect  gas 
with  constant  specific  heats  and  the  perturbation  method 
developed  by  Professor  Hui  is  employed  to  calculate  the 
resulting  unsteady  flows  over  the  upper  and  lower  surfaces 
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of  a  two-dimensional  flat  plate.   Finally,  the  strip  theory 

approximation  is  utilized  to  combine  the  effects  of  these 

flows  for  the  case  of  a  three-dimensional  wing,  provided  that 

the  bow  shock  is  attached  and  therefore  the  flows  are  independent, 

The  unsteady  flows  over  the  lower  and  upper  sides  of  a 

two-dimensional  flat  plate  are  studied  in  [Ref .  6,7]  and 

[Ref.  8]  respectively.   Nevertheless,  we  will  indicate,  in 

the  following  first  two  sections,  III.A.l  and  III. A. 2,  the 

way  in  which  these  flow  problems  are  formulated  and  solved. 

We  will  also  give  the  solutions  for  the  complete  set  of  flow 

quantities  in  the  disturbed  regions.   These  flow  quantities 

will  be  used  in  Section  IV  where  the  effects  of  upstream 

unsteadiness  in  the  flow  are  considered. 

1 .   Two-Dimensional  Oscillating  Flat  Plate — Compression 
Side 

Instead  of  a  flat  plate  the  equivalent  flow  problem 

of  a  two-dimensional  wedge  is  considered.   This  problem  is 

formulated  and  solved  in  [Ref.  6,7]  with  the  ultimate  goal 

of  studying  the  stability  of  wedges/caret  wings.   In  what 

follows  in  this  and  the  next  subsection, 

a)  The  major  steps  in  the  method  of  solution  are 
indicated. 

b)  A  generalized  approach  that  permits  the  formulation  of 
the  fluid-shock-fluid  boundary  conditions  is  adopted.   This 
approach  is  described  in  Appendix  A  and  the  formulation  of 
the  boundary  conditions  for  the  two  cases  is  given  in  Appen- 
dices B  and  C. 
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c)  Some  of  the  results  and  discussions  contained  in 
[Ref.  6,7,8],  pertaining  to  the  stability  of  wedges/caret 

wings  and  flat  plates,  are  not  included  since  they  are  not 
directly  related  to  the  subject  matter  of  this  thesis. 

d)  On  the  other  hand  the  solutions  for  the  complete  set 
of  the  flow  quantities  in  the  disturbed  regions,  which  are 
not  included  in  the  above  references,  are  given.   Much  of 
the  mathematical  detail  in  obtaining  these  solutions  is 
omitted.   These  flow  quantities  will  be  used,  as  mentioned 
before,  in  Section  IV  where  the  effects  of  upstream  unsteadi- 
ness in  the  flow  are  considered. 

e)  The  same  symbols  as  those  used  in  the  references  will, 
in  general,  be  employed.   Changes  will  be  limited-  to  those 
necessary  for  clarification  purposes  or  generalization  of 
approach. 

a.   Problem  Formulation 

Consider  a  two-dimensional  wedge  of  length  1 ,  at 
design  condition  (zero  mean  angle  of  attack) ,  in  a  supersonic/ 
hypersonic,  uniform,  steady  flow  of  an  inviscid  perfect  gas 
with  constant  specific  heats  (Figure  2.1a) .   Assume  that  the 
wedge  is  performing  a  low  amplitude  and  frequency  harmonic 
oscillation  in  pitch  with  given  circular  frequency  oj,  about 
an  axis  perpendicular  to  the  plane  of  the  paper,  through  the 
point  C  shown.   Let  a  system  of  cartesian  coordinates  Oxy 
be  attached  to  the  wedge  so  that  0  is  at  its  apex  and  axis 
Ox  is  along  the  mean  position  of  the  upper  surface.   The  bow 
shock  is  assumed  to  be  attached  to  the  body  and  the  flow 
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wedge  surface  in 
steady  flow 

unsteady  fbw 

shock  surface  in 

steady  reference  flow 
unsteady  flow 


Figure  2.1a.   Oscillating  wedge  in  uniform  steady  flow 


flat    plate 


shock 


Figure  2.1b.   Oscillating  flat  plate— Compression  side 
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quantities  on  the  upper  surface  of  the  wedge  are  to  be 
found.   For  all  practical  purposes  the  solution  to  this 
problem  will  give  the  flow  quantities  on  the  lower  side  of  a 
two-dimensional  flat  plate  of  length  I    sec  9  at  an  angle  of 
attack  9  (Figure  2.1b) . 

b.   Method  of  Solution 

The  unsteady  flow  over  the  upper  surface  of  the 
wedge  will  be  found  by  perturbing  the  steady  shock  flow  con- 
sidered in  Section  II. A. 4. a. 

Denote  by  U  ,  p  ,  p   the  velocity,  pressure  and 


00       00 


density  in  region  A.   Denote  by  u  ,  p  ,  p   the  velocity, 

pressure  and  density  of  the  steady  reference  flow  in  region 

B.   Non-dimensional  lengths  and  time  are  introduced,  defined 

by 

—  —  u 

x  =  ~,     y   =  -2L   and   t  =  -2.  t       (2.4) 

II  I 

Assume  that,  as  a  result  of  the  oscillation  of 
the  wedge,  the  perturbed  flow  quantities  in  region  B  are 
given  by 


u=u+eu+...  (2. 5a) 

o 


v   =   ev  +  .  .  .  (2 .5b) 


p   =  p   +  ep  +  . . .  (2.5c) 


p   =   p   +  ep  +  .  .  .  (2.5d) 
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where  z   describes  the  deviation  of  the  unsteady  flow  from 
the  reference  steady  flow. 

Substitute  quantities  (2.5)  in  the  governing 
equations  of  motion  (2.1),  non-dimensionalize  the  indepen- 
dent variables,  by  using  relations  (2.4)  and  get  the  following 
perturbation  equations 


ut  +  ux  =   -  ;  „   p„  (2.6a) 


P    u 
o   o 

1 

py 

P    u 
o   o 

vfc  +  vx   =   -  :~r-  p..  (2.6b) 

p.  +  p    =   a  (p.  +  p  )  (2.6c) 

rt   *x      o   t   x 

p 

p.  +  p   +  — (u  +v  )   =   0  (2.6d) 

t   Kx   u   x   y 
o      2 

where  subscripts  denote  partial  differentiation  and  a   is  the 
speed  of  sound  in  the  reference  steady  flow.   Assume  that 
the  perturbation  quantities  u,  v,  p  and  p  have  the  form 


u  =  u  elktU(x,y)  (2.7a) 


v  =   u  elktV(x,y)  (2.7b) 


p   =   p  yM  elktP(x,y)  (2.7c) 


p   =   p  M  elktR(x,y)  (2.7d) 


where  U,  V,  P  and  R  are  unknown  quantities  to  be  found,  M 
is  the  Mach  number  in  the  reference  steady  flow  and  k  is 
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the  so-called  reduced  frequency  of  oscillation  defined  by 
k  =  a)£/u  .   k  is  assumed  to  be  small  and  the  time  indepen- 
dent quantities  U,  V,  P  and  R  are  expressed  as  power  series 
in  (ik)  of  the  form 


U   =   U(  }  +  (ik)U(1)  +  ...  (2.8a) 


V   =  V(0)  +  (ik)V(1)  +  ...  (2.8b) 


P   =   p'0)  +  (ik)P(1)  +  ...  (2.8c) 


R   =   R(0)  +  (ik)R(1)  +  . . .  (2.8d) 


Expressions  (2.8)  are  substituted  in  (2.7)  and  the  resulting 
expressions  in  the  perturbation  equations  (2.6).   By  equating 
the  same  order  terms,  in  each  of  these  equations,  a  sequence 
of  systems  of  partial  differential  equations  is  formed. 
Only  the  systems  of  zeroth  and  first-order  equations  are  of 
interest  in  stability  analysis.   The  zeroth-order  equations 
are 


0<o>   .     lp(0) 
x  M    x 


(2.9a) 


o 


V 


(0)   =   -  -i-P(0)  (2.9b) 

x         My 
o   J 

P(0)   =   R(0)  (2.9c) 

X  X 


u<o>  +  v(0)  +  M  R(o>  _  (2-9d) 

x      y      ox 
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The   first-order   equations    are 


0(i>    +ip'l>      .     _<«  C2.10a, 

O 

v'1'    .ip'1'      =     -V<0»  (2.10b) 

pd)  .  R(i)     .    r(o)   .  p(o)  (2-10c) 

XX 


0(1)    +   V(1)    +   M   R(1)       =      -M   R(0)  <2.10d) 

x  y  o    x  o 


Next  we  consider  the  boundary  conditions  applicable 
to  the  problem. 

Along  the  surface,  the  condition  to  be  satisfied 
is  the  flow  tangency  condition  given  by  equation  (2.2). 
The  equation  of  the  surface  for  a  stationary  wedge  is  given 

by  S  =  y  =  0  and  for  an  oscillating  wedge  is  given  by 

*  ikt 
S(x,y,t)  =  y  +£ (h  cos  9  -  x)  =  0  where  e  =  ee    .   Equation 

(2.2)  then,  gives  with  V  =  {u,v}  =  u  {1  +£U,eV}, 

V(x,y)   =   1  +  (ik)(x-h  cos  9)  at  y  =  e(x-hcose) 

Expanding  V(x,y)    about  y  =   0    and   neglecting   higher-order 
terms   we   get   the   linearized   condition 

V(x,0)      =      1   +    (ik)fc-h  cos  8)      at     y  =   0 

Use   of  equation    (2.8b)    gives 

V(0)    =1  At     y      =      0  (2.9e) 
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V(1)   =   x-hcose    At  y   =   0  (2.10e) 


Across  the  shock,  the  conditions  to  be  satisfied,  are  the 
conditions  given  by  (2.3).   The  equation  of  the  shock  in 
the  steady  reference  flow  is  given  by  G   =  -y  +  x  tan  <j>  =  0 . 
Let  the  equation  of  the  shock  in  the  unsteady  flow  be  given 
by  G  =  -y  +  x  tan  <f>  +  sQ(x)  =  0,  where  Q(x)  is  an  unknown 

function  to  be  determined  as  part  of  the  solution  and 

A  ikt 
£  =  ee  "'  .   To  find  the  boundary  conditions  across  the  shock 

we  substitute  expressions  (2.7)  in  (2.5)  and  the  resulting 

expressions  in  equations  (2.3).   The  boundary  conditions, 

after  linearization,  are  given  by 


V  =  AQ'  +  (ik)BQ  AT  y  =  x  tan  <j>          (2.11a) 

P  =  CQ'  +  (ik)DQ  (2.11b) 

U  =  EQ1  +  (ik)FQ  (2.11c) 

R  =  GQ'  +  (ik)JQ  (2. lid) 

where  Q'  =  dQ/dx  and  the  values  for  the  constants  A  through 
J,  which  depend  on  the  reference  steady  flow,  are  given  in 
Appendix  B.   The  derivation  of  these  relations  is  lengthy 
and  tedious,  even  for  the  case  considered  here,  where  there 
are  no  upstream  disturbances.   In  Appendix  A  equations  (2.3) 
are  put  in  an  alternate  form.   This  form  permits  a  much  easier 

solution  of  the  equations  and  is  repeatedly  used  throughout 

o 
this  thesis. 
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Expressing  the  time  independent  quantities,  U, 
V,  P  and  R  by  their  expressions  as  given  by  equations  (2.8) 
we  get  the  following  boundary  conditions  for  the  zeroth  and 
first-order  systems  respectively. 


At  y  =  x  tan  <f>   V(0)   =   AQ '  (0) 


P(0»   =   CQ 


U(0»   =   EQ 


R(0)  =    aQ 


at  y  =  x  tan  <p        V 


(1) 


U(l) 


=  AQ 


P(1)   =   CQ 


=   EQ 


R(1)   =   ^ 


(0) 
(0) 
(0) 


(1>  +  BQ(0) 


(1>  +  DQ(0) 


(1)  +pQ(0) 


(1)  +HQ(0) 


(2.9f) 

(2.9g) 

(2.9h) 

(2.9i) 

(2.10f) 

(2.10g) 

(2.10h) 

(2.10i) 


c.   Complete  Solution 

Two  boundary  value  problems  have  been  set  up. 

The  zeroth-order  equations  (2.9a-2.9d)  and  the 

zeroth-order  boundary  conditions  (2.9e-2.9i)  constitute  the 

first  boundary  value  problem.   This  problem,  which  will  be 

solved  first,  corresponds  to  the   problem  of  steady  flow  past 

a  wedge  and  its  solution  should  give  the  flow*  quantities  .... 

behind  the  shock  for  a  stationary  wedge.   This  result  is 

shown  in  Appendix  B. 
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The    first-order  equations    (2 . 10a-2 .10d)    and   the 
first-order  boundary   conditions    (2.10e-2 .10i)    constitute 
the    second  boundary  value   problem. 

In  both  problems    the   equations    and  boundary   con- 
ditions   are   linear   and   therefore   suggest   for   the  unknowns 
U(l),    V(l),    P(l)    and   R(l),     (i    =   0,1),    solutions    that  are 
linear   combinations   of    the   non-dimensional   spatial   coordinates 
x   and  y. 

In   view  of    the   above   reasoning,    we   assume    for   the 
first  problem  a   solution  of   the   form 

U(0)  =  u<°>x   +  u<°>y  +  u<°>  (2.12a) 

V<0>  =  v<°>x   *   v<0)y  +  v<0)  (2.12b) 

P(0)  =  p{°>x  +  p<0)y  +  p<°>  (2.12c) 

R(0)  =  r<0)x   +   r<°>y  +  r<0)  (2.12d) 

Q<°>  =  q<°>x   +  q<°>  <2.12e) 


Substitution   of   these   values    into    the    zeroth-order   equations 
and  boundary  conditions    (2.9)    gives    the    following   solution 


U(0)      =      u<0)       =      l/l  (2.13a) 


V(0)      =     v^0)      =      1  (2.13b) 

P«»       =      p^°)      =      c/l  (2.13c) 
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R(0)  m      r^0)   _   Q/A  (2.13d) 

Q(0)   =   x/A  +  q^0)   =   x/A  -  h  cos  9         (2.13e) 

where  the  coefficient  q^   was  found  from  the  additional  con- 
dition  that  the  bow  shock  is  attached  to  the  leading  edge 
of  the  body.   Recalling  that  the  equations  of  the  surface 
and  the  shock  are  given  by  S(x,y,t)  =  -y+e (x-hcos6)  =  0 
and  G(x,y,t)  =  -y+xtan^+eQ (x)  =  0  we  get  the  above  result 
by  letting  x  =  0  in  the  equation  e(x-hcos6)  =  xtan<J>+£  (Q 
+  (ik)Q(1))  . 

Similarly  to  solve  the  second  boundary  value 
problem  we  assume  a  solution  of  the  form 

U(1)  =  uf»  x  +  uf'y  ♦  u*1'  (2.14a) 

V(1)  =  v^x  ♦  v<X)y  +  v<X)  (2.14b) 

P(1»  =  pf'x  +  p<X)y  +  e'1'  (2.14e) 

R(1)  =  r™x  *   r«X,y  +  r<X)  (2.14d) 

Q(1»  =  q^'x2  +  q^'x  +  q*1'  (2.14.) 


and  substitute  these  values  to  the  first-order  equations  and 
boundary  conditions  (2.10)  .  The  resulting  solution  is  given 
below 


uj1)   =   -[Cq^/M  +  (E  +D/Mq)/A  +  2tan<j>] 
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ui1}       =       [q.(1)  (E+C/M    )+(E+F+D/M    ) /A] /tan  <J>  +2 
«  J.  o  o 

u^       =      h  cos  9  [E(B-1)/A -F] 


v<X)      =      1,  v^      =      -hcos 


v^X)       =       (Aq;[1)    +B/A-l)/tan    0 


p|1)       =      Cq^1*    +   D/A  +    2M      tan   <j> 


p'11       =       -    2   M 

r2  o 


p^1)       -      h   cos    6  [C(B-1)/A  -D] 


r^1'       =      Cq|1)    +    (D+C-G)/A  +    2  M      tan   <j> 


r^1}       =       [(G  -C)q;[1)   +  (G  +H  -C  -D)/A]/tan   <J>-  2M 
r.j1}       =      h    cos    9  [G(B-1)/A -H] 

q{1}       =       [k2(A-B-C+E/M    ) /A  -M     tan  4>[D/A+2  M    tan<j>)  ]  /2  (k2A 

+  CM     tand>) 

o  T 

2  2         2 

k^      =      MV(M-l) 
o        o 


2.   Two-Dimensional  Oscillating  Flat  Plate — Expansion 
Side 

This  problem  is  formulated  and  solved  in  [Ref .  8] . 
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a.  Problem  Formulation 

Consider  a  two-dimensional  flat  plate  of  length 
I  in  a  supersonic/hypersonic  uniform  steady  flow  of  an 
inviscid  perfect  gas  with  constant  specific  heats.   Assume 
that  the  plate  is  performing  a  low  amplitude  and  frequency 
harmonic  oscillation  about  its  leading  edge.   Let  a  system 
of  Cartesian  coordinates  Oxy  be  attached  to  the  body  so  that 
0  coincides  with  the  leading  edge  of  the  plate  and  axis  Ox 
is  along  its  mean  position  (Figure  2.2a).   The  bow  shock  is 
assumed  to  be  attached  to  the  leading  edge  and  the  flow 
quantities  on  the  upper  surface  of  the  plate  (in  Region  C) 
are  required. 

b.  Method  of  Solution 

The  unsteady  flow,  over  the  upper  surface  of 
the  oscillating  plate,  will  be  found  by  perturbing  the 
steady  Prandtl-Meyer  flow,  considered  in  Section  II.A.4.b. 
Denote  by  U  ,  p   and  p   the  velocity,  pressure  and  density 
in  Region  A.   Denote  by  u.  ,  p..,  p,  the  velocity,  pressure 
and  density  of  the  reference  steady  flow  in  region  C.   The 
solution  procedure  that  follows  is  similar  to  the  one  pre- 
sented in  Section  III.A.l.b  and  most  of  the  assumptions  and 
results  given  there,  apply  to  this  section  too,  provided 

that  u  ,  p  ,  p  ,  M   and  a   (the  reference  steady  flow  quanti- 
o   o   o   o      o 

ties  over  the  compression  side  of  the  flat  plate)  are  replaced 
by  u.  ,  p1  ,  p..  ,  M..  and  a,  (the  reference  steady  flow  quantities 
over  the  expansion  side  of  the  flat  plate). 
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U«a,    M 


Front     in 
-.  .  unsteady  flow 
reference    steady  flow 


Figure  2.2a.   Oscillating  flat  plate — Expansion  side 
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Figure  2.2b.   Oscillating  flat  plate— Polar  coordinates 
for  P r and tl -Meyer  flow 
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Non-dimensional  time  and  lengths  are  introduced, 
defined  by  (2.4).   As  a  result  of  the  oscillation  of  the 
plate  assume  that  the  perturbed  flow  quantities  in  region  C 
are  given  by  (2.5).   The  resulting  perturbation  equations  are 
given  by  (2.6).   Assume  that  the  perturbation  quantities  are 
given  by  (2.7)  and  let  the  time  independent  quantities  U,  V, 
P  and  R  be  expressed  as  power  series  in  (ik)  by  (2.8)  .   The 
zeroth  and  first-order  equations,  derived  as  before,  are 
given  by  (2.9a-2.9d)  and  (2 . 10a-2 . lOd)  restated  below. 


U(0»  =   -^P(0)                           (2.16a) 

x  M,   x 

V<0)  =   -  iP<°)                             (2.16b) 

x  M,   y 

P(0)  =   R(0)                               (2.16c) 

X  X 

U(0)  +  V(0)  +  M.  R(0)   =   0                  (2.16d) 

x  v       1   X 


0d)  +  i  p(i)  .  .  o(0)              (2>17a) 

x  M,   x 

V'1'  +  1-p'1'   =   -  v(0)                    (2.17b) 

x  M,   y 

p(D  .  (i)  _   (o,    (0)            (2<17c) 

X  X 

U(1)  +  V(1)  +Ml  R(1)   =   -  Mn  R(0>           (2.17d) 

x  y      1  x         1 


We  now  consider  the  boundary  conditions  applica- 
ble to  the  problem. 
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Let  the  equation  of  the  oscillating  surface  of 
the  plate  be  given  by  S(x,y,t)  =  y  -ex  =  0.    The  flow 
tangency  condition  (2.2),  with  V  =  u,{1+eU,£V},  gives  after 
linearization 


V(0)   =1         At  y   =   0  (2.16e) 


V(1)   =   x         At  y   =   0  (2.17e) 


We  next  consider  the  boundary  conditions  across 
the  surface  that  is  separating  regions  B  anc  C.   We  assume 
that  as  a  result  of  the  small  amplitude  slow  oscillations 
of  the  body  the  separating  surface  is  slightly  deformed  and 
its  equation  is  given  by  G  =  -y  + x  tan  <p   +eQ(x)  =  0  where 
Q(x)  is  an  unknown  function  which  may  be  expanded  as 
Q(x)  =  Q    +  (ik)Q   '+....   We  call  the  flow  expansion 
an  expansion  front  or,  simply,  front  and  assume 
that  upstream  of  it  the  Prandtl-Meyer  flow  is  not  disturbed 
while,  along  the  front,  the  unsteady  flow  matches  the  steady 
Prandtl  Meyer  flow  continuously.   The  assumptions  made  are 
completely  analogous  to  the  assumptions  made  in  the  case  of 
a  finite  compression  shock  discontinuity  and  the  Rankine- 
Hugoniot  conditions  (2.3)  may  be  used  to  give  the  boundary 
conditions  across  the  expansion  front.   We  note  that,  since 


The  difference  in  the  form  of  this  equation  and  the  one 
considered  in  the  previous  section  is  due  to  our  assumption 
that  in  this  case  the  flat  plate  is  oscillating  about  its 
LE  and  thus  h  =  0. 
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we  are  dealing  with  an  infinitely  weak  discontinuity,  the 
flow  across  the  front  is  isentropic.   The  procedure  is  given  in 
Appendix  C  and  the  resulting  boundary  conditions,  after  lineari- 
zation and"  use  of  equations  (2.8),  are  given  below. 


V(0)   =   A'Q'  (0)  At  y   =   xtan<J>     (2.16f) 


(0)   =   C'Q'(0)  (2.16g) 


U(0)   =   E'Q'(0)  (2.16h) 


R(0)   =  G'Q'(0)  C2.161) 


V(1)   =   A'Q1  (1)  +  B'Q(0>   At  y=xtan  (j>      (2.17f) 
p(D   =   C,Q,  (1)  +  d.q(°)  (2.17g) 

U(1)   =   E'Q'  (1)  +  F'Q(0)  (2.17h) 

R(1)   =   G'Q'(1)  +  J'Q(0)  (2.17i) 

where  the  coefficients  A'  through  J'  depend  on  the  reference 
steady  flow  and  are  given  in  Appendix  c. 
c.   Complete  Solution 

Two  boundary  value  problems  have,  again,  been 
set  up  and  will  be  solved  successively.   The  zeroth-order 
equations  and  boundary  conditions  (2.16)  constitute  the 
first  problem,  which  corresponds  to  a  steady 
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Prandtl-Meyer  flow  problem.   To  solve  it  assume,  as  before, 
that  the  unknowns  have  the  form  given  by  (2.12),  substitute 


in  (2.16)  and  get  with  <  =  M ./  VM2  -1 

P<0)   =  p<0)   =  <  (2.18a) 

R<0)   =  r<0)   =  k  (2.18b) 

U(0)   =  u<°»   =  -K/Ml  (2.18c) 

V(0)   =  v<°>   =   1  (2.18d) 

Q«0)   =  I+i  K4  x  +  q<°>  (2.18e, 

where  the  coefficient  q^    =  0  since,  in  this  case,  the 
plate  is  oscillating  about  its  LE . 

The  second  boundary  value  problem  consists  of 
the  first-order  equations  and  boundary  conditions  given  by 
(2.17).   In  this  case  the  unknowns  are  expressed  by  (2.14) 
and  the  resulting  solution  is 

(i)        VMi"2) 

Pu;   =   L      L             x  -  2Mxy  (2.19a) 
(  v/m^-1)  3 

R(1)   =   P(1)  (2.19b) 


(2.19c) 


0(1) 

= 

1 

3M^  +1 

M^  -1 

(    v/m^-1)3 

0 
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(1)         2M*-1 
VK±)      =   x x     y  (2.19d) 

>/m^  -1 

3.   Three -Dimensional  Oscillating  Wings  of  Arbitrary 
Planform  Shape 

The  problem  is  formulated  and  solved  in  [Ref.  5] 
and  only  the  major  steps  in  the  method  of  solution  are  given 
here. 

a.   Problem  Formulation 

Consider  an  oscillating  wing  of  arbitrary  plan- 
form  shape  at  an  angle  of  attack  in  a  steady,  uniform,  super- 
sonic/hypersonic flow  (Figure  2.3).   Assume  that  the  oscilla- 
tions are  periodic  with  small  amplitude  and  frequency  and 
that  the  bow  shock  is  attached  to  the  wing.   We  let  the 
pressure,  density,  velocity  and  Mach  number  of  the  approach- 
ing flow  be  given  by  p^,  p  ,  U^  and  M^ .   We  also  denote  the 
total  area  and  the  root  chord  of  the  wing  by  S  and  I   respec- 
tively and  the  distance  of  the  pivot  position  from  the  leading 
edge  by  x  .   We  assume  that  the  pitching  motion  of  the  wing 
is  described  by 

9(t)   =   6eiwt 

where  9~  and  u  are  the  amplitude  and  circular  frequency  of 
oscillation  and  t  is  the  non-dimensional  time. 

We  define  the  reduced  frequency  of  oscillation 
by 


k   =   oo£/U 

'    00 
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x/l  =  f(y/b) 


x/l  =  g(y/b) 


Figure    2.3.      Three-dimensional  wing 
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b.   Method  of  Solution 

The  in-pitch  stability  derivatives  of  the  wing 
are  required.   The  pitching  moment  coefficient  C  ,  the 
stiffness  derivative  -C   and  the  damping  in-pitch  deriva- 

me 

tive   -C        are   defined   by 

Cm      =      2M/pooU^S      =       (IAS)   J  J(x-xc)C    (x,y,t)dS 


0(t) [(-C      )   +  (ik) (-C      ) ]  (2.20) 

m0  m,? 


where  M  is    the  moment   about   the  pivot   axis    and  C      is    the 

P 

pressure  coefficient. 

The  pressure  coefficient  is  defined  as  usual  by 

2(P-Pj 
C 


P      p  u2 

r  CD      00 

and  using  the  two-dimensional  flow  assumption  can  also  be  written  in  the  form 

x 

C        =       (C)       +    9(t)[A+ik(Bi-c4)] 
p  p  x.  Jc 

r  v  o 

=       (C)       +    9(t)  [A+i^(Bx  -CxJ]  (2.21) 

P    o  Uco  c 

where  A,  B  and  C  are  dimensionless  constants  which  are 
functions  of  geometry  and  steady  flow  quantities  and  (C  ) 
is  the  mean  pressure  coefficient  corresponding  to  9 (t)  =  0. 
This  coefficient  has  no  contribution  to  the  stability 
derivatives  and  will  be  neglected. 

Our  goal  is  to  calculate  the  coefficients  A,  B, 
C  appearing  in  relation  (2.21)  and  use  the  pressure  coefficient 
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in  relation  (2.20)  to  obtain  the  stability  derivatives. 
It  is  clear  that,  for  a  two-dimensional  flat  plate  with 
attached  shock  wave,  the  coefficients  will  be  the  sum  of 
independently  obtained  coefficients  over  the  compression 
and  expansion  sides  of  the  plate  respectively.   Thus  the 
coefficients  will  have  the  form 


A  =   A.  -  A  ,    B   =   B„  -  B  ,   C   =   C„-C      (2.22) 
£     u  I  u  I  u 


where  the  subscripts  I   and  u  denote  the  lower  and  upper 
surface  respectively. 

The  results  from  Section  III.A.l  and  [Ref.  6,7] 
are  used  to  obtain  coefficients  A,,  B,  and  C,  of  the  form 


A^   =   AoC/A,    B£   =   yo(2G-I),    C%      =   UQI  (2.23a) 


where  X      and  u   are  coefficients  introduced  to  account  for 
o      o 

differences  in  notation  given  by 

O   oo     oo 

and  G,  I  are  quantities  defined  in  [Ref.  6,7]. 

Similarly  the  results  from  Section  III. A. 2  and 
[Ref.  8]  are  used  to  obtain  coefficients  A  ,  Bu  and  Cu  of 
the  form 


Au   =   -A1M1/(m2-1)0-5,     Cu   =   Vl/Xl  (2.23b; 

Bu   =   uiM1(mJ-2)/(M^-1)1*5  , 

78 


where  A.  ,  y,  are  correction  coefficients  to  match  present 
notation  of  the  form 

2  pl   ul  2 

00       00 

c .   Solution 

Substituting  relations  (2.23)  in  relation  (2.22) 
and  putting  the  resulting  expressions  in  relation  (2.21) 
we  obtain  the  pressure  coefficient  for  a  two-dimensional 
flat  plate.   We  shall  employ  the  strip  theory  to  solve  the 
three-dimensional  wing  problem  in  hand,  which  means  we  will 
assume  that  the  flow,  at  each  point  of  the  wing,  is  two- 
dimensional  locally.   This  assumption  permits  the  use  of  the 
two-dimensional  flat  plate  pressure  coefficient  for  the  case 
of  the  three-dimensional  flat  wing. 

Using  relations  (2.21)  and  (2.20)  we  get 


-C     =   A(I.  -  x  /£)  (2.24a) 

mQ        1    c 

-C  =       [BI,+  (C-B)I,]-[(B+C)I,  +  (C  -B)I.]x   /l+   CU/i)2 

m  z  **  ±  j      c  c 


(2.24b) 


where 


I.   =  k   /  (g2-f2)dn,     I,   =   2k/3   /  (g3-f3)dn   (2.25) 
1        0  ^  0 


1  1    2   2 

I.  =  2k.   /  f(g-f)dn,   I.  =  k  /  f(g  -f  ) dn 
J       x  0  0 

k  =  ib/S,  Q  n  =  y/b 
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For  delta  wings  with  power  law  leading  edges,  i.e.,  wings 
with  f (n)  =  n    and  g(n)  =1,  equations  (2.25)  become 


I1      =   (n+l)/(n+2),      I2   =   (n+l)/(n+3)        (2.25a) 
I3   =   n/(n+2)  ,         l4   =   n/(n+3) 


The  minimum  value  of  C   obtained  for  pivot  position 

m9 

X  /£   =   [B+(2n+l)C]/2C(n+2) 


is 


(-C.)  .    =  -£CB*(2!£»C  l2  +  4(B+nC)/(n+3) 

in  mm       c     n+^ 


By  setting  (-C   )  •   =  0  a  stability  boundary  for  power  law 

nu  min  ■*        ■*     r 

delta  wings  may  be  obtained,  which  is  practically  independent 


of  the  power  n,  as  shown  in  (Figure  2.4) .   Plots  of  C   and 

mQ 

C   vs  a  for  several  values  of  pivot  axis  position  and  powers 
m9 

n  are  given  in  (Figure  2.5).   In  (Figures  2.6,  2.7)  the 
stability  derivatives  vs  the  pivot  position  are  plotted  for 
several  values  of  power  n  and  angles  of  attack  10°  and  20°. 
Comparisons  of  results  obtained  by  the  present  theory  with 
results  obtained  by  other  theories  and  related  discussions 
are  included  in  [Ref.  5]  and  will  not  be  repeated  here. 
Nevertheless  in  Section  III.C  we  will  compare  the  present 
results  with  potential  flow  theory  results. 
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Figure  2.4. 

Stability  boundary 
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Figure  2.5.   Stability  derivatives  vs  AOA 
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Figure  2.6.   Stiffness  derivative  vs  pivot  position 
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Figure  2.7.   Damping  derivative  vs  pivot  position 
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B.   LINEARIZED  THEORY 

It  was  pointed  out  in  Section  II  that  several  methods 
are  used  in  the  linearized  theory  of  supersonic  aerodynamics. 
A  simple  application  of  one  of  these  methods,  the  so-called 
method  of  fundamental  solutions,  will  be  described  in  this 
subsection.   For  this  application  a  distribution  of  pulsating 
sources  over  a  "simple"  planform  will  be  considered. 

Before  dealing  with  this  application,  however,  the  notion 
of  a  disturbance  propagation  in  supersonic  flow  will  be  re- 
viewed and  the  fundamental  solution  for  a  moving  pulsating 
source  will  be  introduced. 

1.   Propagation  of  Disturbances 

Supersonic  flow  is  dominated  by  the  fact  that  dis- 
turbances travel  with  finite  velocities,  namely,  the  speed  of 
sound.   In  formulating  the  linearized  potential  equation 
(1.28)  deviations  of  the  speed  of  sound  from  its  free  stream 
value  a  were  neglected  and  also  the  perturbation  velocities 
u1 ,  v1 ,  w1  were  taken  very  small  compared  to  the  free  stream 
velocity  U^.   As  a  result  any  disturbance  from  a  source 
located  at  point  (x,y,z) ,  in  a  coordinate  system  fixed  to 
a  body  in  the  flow,  can  be  felt  only  inside  or  on  the  surface 
of  a  right  cirular  cone  whose  axis  points  downstream  from 
the  source  (Figure  2.8) .   An  observer  moving  with  the  fluid 
sees  a  pulse  emitted  at  t  =  0 ,  expanding  on  a  spherical  sur- 
face with  instantaneous  radius  a  t  and  center  moving  down- 
stream  with  velocity  U  .   The  positions  of  these  expanding, 
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Upstream   zone    of    influence 


Downstream   zone    of  influence 


Figure  2.8.   Upstream  and  downstream  zones  of  influence. 
Disturbance  located  at  (£,n*?)  at  time  t  =  0 

moving  spheres  form  an  envelope  which  is  a  cone  of  semivertex 
angle 


.  -i  a»t 
p   =  ain    _ 

00 


.  -1    1 

sln  ST" 

M 

00 


=      tan 


-1 


\/m2-i 


where  u  is  the  Mach  angle. 

This  cone  is  known  as  the  Mach  cone  or  downstream 
zone  of  influence  of  the  point  (x,y,z)  and  its  equation  is 
given  by 


(C-x)2  -  (m2-1)  [(n-y)2  +  U-z)2]   =  0       (2.26) 


On  the  other  hand  the  point  (x,y,z)  can  be  influenced 
by  sources  whose  locus  is  evidently  a  similar  cone  directed 
forward  from  (x,y,z) .   This  cone  is  known  as  the  forecone 
or  upstream  zone  of  influence  and  has  the  same  equation  except 
that  for  this  case  x  >  E, . 
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Assume  now  that  a  steady  supersonic  flow  has  been 
established  past  a  stationary  three  dimensional  wing  lying 
very  close  to  the  xy-plane.   Regarding  each  point  of  the  wing 
as  a  disturbance  source  we  can  see  that  the  downstream  zone 
of  influence  for  the  entire  wing  is  bounded  by  the  envelope 
of  the  Mach  cones  emanating  from  the  leading  edge.   If  the 
leading  edge  is  straight  (or  the  wing  is  two-dimensional) 
the  envelope  reduces  to  the  so-called  Mach  wedges. 

To  calculate  the  fluid  motion  at  any  point  (x,y,z) 
we  need  to  consider  only  the  contribution  from  the  disturbance 
sources  that  belong  to  the  region  of  the  xy-plane  intercepted 
by  the  forecone  from  point  (x,y,z).   This  area  of  influence 
forms  a  hyperbola  and,  with  sources  assumed  to  lie  on  the 
xy-plane,  is  found  from  the  equation  of  the  forecone  (2.26) 
by  setting  <;  =  0 .   We  thus  get 


HX  2  =     y   ±      ^    (x-02/(Mf-l)-z2  (2.27) 


A  generalized  supersonic  planform  with  leading  edge 
AA'C'C,  trailing  edge  DD'F'F  and  streamwise  tips  AF  and  CD 
is  considered  next  (Figure  2.9) .   For  each  point  (x,y,0  )  the 
area  of  the  sources  that  influence  the  point  reduces  to  a 
region  bounded  by  two  straight  lines  upstream  of  the  point. 
These  lines  found  from  the  forecone  equation  (2.26)  by  setting 
z  =  0  and  z,   =  0  are  given  by 


£-x  =  ±  vM  -1  (y-n)     C  :.  x 
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Limit   of 
undisturbed   flow 


Figure  2.9.   Generalized  supersonic  planform 

and  are  shown  as  dotted  lines  making  an  angle  y  with  the 
x-axis  (Figure  2.9). 

The  portions  A'BC  and  D'EF'  of  the  leading  and 
trailing  edges  are  called  supersonic  since  the  velocity  normal 
to  these  edges  is  greater  than  the  speed  of  sound.   Similarly 
the  remaining  portions  of  the  leading  and  trailing  edges  are 
called  subsonic  since  the  normal  component  of  velocity  is 
less  than  the  speed  of  sound.   Along  the  supersonic  portions 
of  the  leading  and  trailing  edges  there  is  no  communication 
between  the  upper  and  lower  surfaces  of  the  wing  and  the  flow 
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over  the  top  or  bottom  of  the  planform  can  be  calculated 
without  reference  to  the  shape  of  the  opposite  side.   On 
the  contrary  the  upper  and  lower  surfaces  of  the  wing  are 
not  independent  along  subsonic  portions  of  the  edges . 

Points  1,  2   and  3  shown  (Figure  2.9)  are  selected 
to  illustrate  different  sorts  of  upstream  influence  regions. 
The  difficulty  in  solving  the  linearized  flow  problem,  i.e., 
finding  the  fluid  motion  at  these  points,  increases  as  we 
move  from  point  1  to  point  3. 

There  exists  no  universal  method  of  approach  in  solving 
the  linearized  problem  for  different  sorts  of  influence  re- 
gions.  Thus  each  planform  shape  calls  for  a  different  method 
of  approach. 

It  is  the  simplest  case  of  a  planform  with  purely 
supersonic  leading  and  trailing  edges,  the  so-called  simple 
planform,  that  will  be  considered  later  in  this  section. 
2.   Fundamental  Solution  of  a  Moving  Source 

For  Uot  =  0  the  linearized  potential  equation  (1.28) 
reduces  to 

y    +   y        +   «f)    = *  w 

xx    yy     zz         2   tt 

00 

This  is  the  classical  wave  equation  for  the  propagation  of 
sound  in  a  still  medium  and  its  fundamental  solution  is  given 
by 

Y(x,y,z,t)   =   |F(t-~)  (2.28) 

00 
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where  affl  is  the  speed  of  sound  of  the  undisturbed  flow  and 

I    2    2    2~ 
r  is  the  radial  distance  from  the  origin,  i.e.,  r  =  vx  +  y  +  z 

A  solution  of  the  linearized  potential  equation  can 
thus  be  readily  found  if  we  can  transform  it  into  the  classi- 
cal wave  equation.   To  achieve  this  the  following  transfor- 
mation, known  as  the  Lorentz  transformation,  is  employed. 

X     =  —  U   x,    n=yy,    "c  =  yz 


y 

A 

-M2 

00 

X, 

ya 

OQ 

A- 

M2 

00 

Moo 

T   =   — —  (t  +  ^-T-  x) 

a  a  (1  -Ml 

00  00 

The  velocity  potential  of  a  sound  source  fixed  with 
respect  to  the  xyz  system  is  then  found  to  be 


ya      vl  -  M  -M   x  +  R 

A  oo  oo  co 

Y(x,y,z,t)       =      £ F[ (t )] 

U R  ,,      ..2, 


a    (1  -M    ) 

/,  2  00  x  oo 

vl  -M 


f~2  2         2         2~  2 

where   R  =     vx     +  (1-M   )  (y     +z    )    and   1-M      is   known   as    the 

00     "*  00 

Prandtl  factor. 

To  find  the  constants  a  and  y  introduced  by  the 
Lorentz  transformation  we  require  that  the  sound  source  should 
produce  constant  sound  flux  independent  of  the  free  stream 
Mach  number  and  we  get 


a   =  a    (1  -M2)  ,         y   =   A  -  m2 

00  00      '  00 


Thus  the  velocity  potential  of  the  moving  source  becomes 
if  in  addition  the  relation  A  =  -  j—   is  used, 
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Y(x,y,z,t) 


-M  x  +R 


4ttR 


F(t  - 


-) 


a  (1  -M  ) 


Expressing  the  above  result  in  a  coordinate  system  moving 
uniformly  with  the  source  located  at  (£,n»0  we  obtain 


Y(x,y,z,t)   =   - 


-i_  P(t— SL) 

4nR  *  [X-      a    ' 

00 


4iRF(t-Tl>     (2-29» 


where 


R  =   /(x-?)2  +  (1  -J?)  [(y  -n)2  +  (z  -  02] 


and 


D      -M(x-S)  -l-R 


a 


a  (1  -M  ) 


By  comparison  of  expressions  (2.28)  and  (2.29)  it  is  seen 
that  the  solution  for  a  moving  source  can  be  obtained  from 
the  solution  for  a  stationary  source  by  replacing  the  ordinary 
distance  r  by  R  in  the  amplitude  and  by  D  in  the  phase.   The 
quantities  R  and  D  are  called  amplitude  and  phase  radii 
respectively. 

A  geometric  interpretation  is  given  in  Figure  2.10. 
At  time  t  a  field  point  Q  and  a  source  0  moving  with  velocity 
U^  in  the  negative  x-direction  are  considered.   For  supersonic 
flow  there  are  two  spherical  waves  passing  through  point  Q 
at  time  t.   These  waves  originated  from  the  source  at  times 
t-  T,  and  t-T  ,  at  which  times  the  source  was  located  at 
positions  P,  and  P2  shown  (Figure  2.10a).   For  subsonic  flow 
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R=0. 


U*t,  =  M„  D 


Figure  2.10a.   Supersonic  source 


U„t  =   M„D 


Figure  2.10b.   Subsonic  source 
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there  is  one  spherical  wave  passing  through  point  Q  that 
originated  at  time  t-T  when  the  source  was  located  at  point 
P  (Figure  2.10b).   From  the  geometry  of  the  figure  we  get 


(x-S-U  T)2  +  (y-n)2  +  (z-O2   =   a2T2         (2.30) 


Solving  this  equation  for  time  T  we  get  for  subsonic  flow 
one  real  positive  solution  and  for  supersonic  flow  two  real 
positive  solutions. 

Physically  we  are  looking  for  the  effect  that  a  dis- 
turbance, originating  at  point  (£,n,£)  at  some  time  t-T, 
will  have  at  some  later  time  t  at  a  point  (x,y,z).   In  this 
sense  the  potential  is  a  retarded  potential.   For  supersonic 
flow  the  disturbance  is  first  felt  at  some  point  (x,y,z) 
after  a  certain  time  T,  has  elapsed.   The  point  (x,y,z) 
penetrates  the  wave  front  of  the  disturbed  region  and  be- 
cause, it  is  moving  at  a  speed  greater  than  that  of  the  wave 
front  it  emerges  from  the  disturbed  region  at  some  later  time 
T~  •   For  subsonic  flow  once  the  point  (x,y,z)  penetrates  the 
wave  front  (after  a  certain  time  T  has  elapsed)  it  will  remain 
in  the  disturbed  region  since  its  speed  is  less  than  that  of 
the  wave  front.   Finally  the  nonexistence  of  positive  roots 
of  equation  (2.30)  should  be  associated  with  an  undisturbed 
region,  i.e.,  with  ¥  =  0. 

In  view  of  above  reasoning  the  source  solution  for 
supersonic  flow  takes  the  form 
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<F(x,y,Z,t)   =   -  JLtFU-T^)  +F(t  -T2)] 


where : 


_     1  M,(x-e)-R 

00  oo       M  —1 


D      1  Mw(x-5)+R 

2     a      a  L    x  2  .   J 

00  °°       M  -1 


R  =   Ax-O2  -  (M2-i)  [(y-n)2  +  (z-c)2] 

For  a  purely  harmonic  time  dependence  the  source 
solution  takes  the  following  final  form 

*<%,y,i,t)   =  -  ^Je-iSr(x-S»cos  A  R  (2.31) 

00 


where : 


co      is    the   frequency  of  oscillation 

oo      =     = is  the  compressible  reduced  frequency  of 

a    (M  -1)       oscillation 


R     =      Mx-5)  2  -  (M2-l)  [  (y-n)  2  +  (z-c)  2] 

3 .   Simple  Planform  Solution 

Consider  a  three  dimensional  wing  performing  a  small 
amplitude  harmonic  oscillation  of  circular  frequency  oj  .   Let 
the  wing  surface  be  very  close  to  the  xy-plane  of  a  Cartesian 
coordinate  system  attached  to  it  and  let  the  body  move  with 
supersonic  velocity  U   in  the  negative  x-direction . 
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The  equations  of  the  upper  and  lower  wing  surfaces 
are  given  by 


Su(x,y,t)   =   z-zu(x,y,t)   =   z  -  hu  (x,y)  ela3t 


S£(x,y,t)   =   z-z£(x,y,t)   =   z  -hJl(x/y)elajt 


and  the  linearized  tangency  conditions  given  by  (1.29)  become 

.  9h 

w(x,y,0    ,t)       =      ela^[icoh     +  U  -j3]  (2.32a) 

u        °°  o  x 

i   1-  3ho 

w(x,y,0~,t)      =      e         [ioih^+U^  -^]  (2.32b) 


Assume  for  simplicity  that  we  are  dealing  with  a 
simple  planform  and  regard  each  point  of  the  wing  as  a  pul- 
sating source.   Recalling  definitions  and  terms  used  in 
Section  III.B.l  this  means  that 

a)  The  leading  and  trailing  edges  of  the  planform  are 
purely  supersonic  and  the  flows  over  the  upper  and  lower 
sides  of  the  wing  are  independent. 

b)  Finding  the  fluid  motion  at  each  point  of  the  surfaces 
involves  the  same  sort  of  upstream  influence  zone,  namely,  the 
sort  indicated  by  point  1  in  Figure  2.9. 

c)  For  a  general  point  (x,y,z)  the  contribution  of  the 
disturbance  sources  that  lie  on  the  hyperbola  with  end  points 
ru  and  n2  given  by  equation  (2.27)  is  to  be  considered. 

In  view  of  the  above  the  complete  solution  for  the 
velocity  potential  can  be  found  by  superimposing  the  effects 
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of   all   sources   that   lie  within   the  upstream  influence    zone 
(Figure    2.11).      Integrating    the   fundamental    solution   for   a 
moving   source,    given  by   relation    (2.31),    over   the   hyperbola 
we  get 


£^x-z    VM2-1 

loot  *  <» 


Y(x,y,z,t)      =      - 


2tt 


?=0 


r       g(£#n)e 


-iuj(x-£) 


cos    rr—  R 

M 

00 


d£  dn 


(2.33) 


where   q(£,n)    is    the   unknown   source    strength. 


forecone 


(x,y,z) 


ZCOtjt  =  zvM^* 


simple    pfonform 


lo) 


X-2<0t|U 


Figure  2.11.   Limits  of  integration 
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To  find  the  source  strength  we  note  that  in  the  pro- 
cedure followed  the  wing  was  represented  by  a  pulsating 
panel  embedded  in  an  otherwise  rigid  xy-plane.   A  mathemati- 
cal boundary  value  problem  has  thus  been  formulated  for  which 
the  normal  velocity  of  the  fluid  on  the  panel  is  prescribed 
by  relations  (2.32)  and  outside  the  panel  it  is  zero,  since 
no  pulsating  sources  exist  there. 

It  can  be  shown  based  on  physical  considerations  that 
the  source  strength  can  be  expressed  in  terms  of  the  normal 
velocity  by  the  realtion 


q(£,n)   =   2w<5,n)   =   2(100*1+1^11) 


Substitution  of  this  value  in  equation  (2.33)  per- 
mits a  straightforward  calculation  of  the  velocity  potential 

Y. 

C.   COMPARISONS  WITH  LINEARIZED  POTENTIAL  FLOW  THEORY 

In  this  section  we  consider  sweptback  tapered  wings  with 
straight  supersonic  leading  and  trailing  edges  and  streamwise 
tips  (Figure  2.12) .   For  this  planform  shape  equations  (2.25) 
become 


Ix   =   £S[(k2-k^)/3  +k2  +1]/S  (2.34a) 

I    =   2£3[(k3-k^)/4  +k^  +3k2/2  +1]/3S  (2.34b) 


I3   =   2£Sk1[(k2-k1)/3  +0.51/S  (2.34c) 
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f(n)=k,n 


M 


(0,0) 


g(n)=K2n  ♦ 


Figure  2.12.   Geometry  of  sweptback  tapered  wings 


2   ,  2 


I        =      £3k1[(k^-kJ)/4  +2k2/3  +0.51/S 


(2.34d) 


S/AB      =      k2"ki+2'  kl      =      xl/l' 


x2/£-l 


Some   special   cases   readily   follow. 
For   a   two-dimensional    flat  plate 


I1      =      0.5,         I2  x/-»#         -3        *4 


=      1/3,         I,  -   IA   -   °'        S/^        =      2 


and  the  stability  derivatives  become 


•C    =  A(0.5  -x  /I) 

mQ  c 


-C     =   B/3-  (B+C)  (x  /JD/2  +C(x  A) 


m 
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For  the  case  of  a  triangular  delta  wing, 
Ix   =   2/3,   I2   =   0.5,   I3   =   1/3,   I4   =   1/4,   SAS   =   1 
and  the  stability  derivatives  are 


-C     =   A(2/3-x/£) 

mQ  c 

-C    =   (B  +C)/4  -  (B  +3C)  (x  A)/3  +  C(x  A) 

me 

For  oscillations  at  small  angle  of  attack  the  coeffi- 
cients A,  B,  C  appearing  in  equation  (2.21)  become 

A   =   C   =   a        B   =   aQ(32-l)/e2         (2.35) 


where  3 


=   7m2-1  and  a   is  the  two-dimensional  lift  curve 
o 


slope  (a   =  4/6) • 

in  this  limiting  case  the  stability  derivatives  given 
above  for  a  two-dimensional  flat  plate  and  a  triangular  delta 
wing  become,  as  expected,  identical  to  well  known  formulas 
based  on  potential  flow  theory  [Ref.  10:  pp.  52,144].   As  a 
result  of  the  strip  theory  approximation  the  stability  deriva- 
tives of  a  rectangular  wing  are  independent  of  its  aspect 
ratio  and  the  stability  boundary  is  always  the  one  for  a 
two-dimensional  wing  [Ref.  10;  Fig.  7.7]. 

Formulas  (2.34)  and  (2.35),  based  on  Hui ■ s  theory, 
were  used  to  calculate  the  stability  derivatives  for  several 
sweptback  wings,  with  straight  leading  and  trailing  edges 
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and  the  results  were  compared  with  results  based  on  potential 
flow  theory,  as  given  in  [Ref .  11,12] .   The  comparisons  are 
shown  in  (Figures  2.13-2.17)  and  the  procedure  followed  to 
read  values  for  the  stability  derivatives  from  [Ref.  11,  12] 
is  described  in  Appendix  E.   Good  agreement  is  generally 
shown . 
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HUI  5  S   THEORY 


4—4-NACA  TN  2699,3196 


A=5.I2,   X=.25,  L  E=3I,  *Q/\  =  0 


Figure  2.13.   Stability  derivatives  vs  Mach  number 
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HUI*S       THEORY 


_^_+NACA    TN     2699,   3196 


LE 


1 0°         20* 


3  0e 


40^ 


50 


»      3» 


60' 


Figure    2.14.      Stability  derivatives   vs    leading  edge 
sweep   angle 
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Figure    2.15.      Stability  derivatives   vs    taper   ratio 
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Figure  2.16.   Stability  derivatives  vs  aspect  ratio 
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Figure  2.17a.   Stiffness  derivative  vs  Mach  number 
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Figure    2.17b.      Damping  derivative  vs  Mach  number 
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IV.   EXTENSION  OF  PROFESSOR  HUI ' S  THEORY — 
UNSTEADY,  NONUNIFORM  UPSTREAM  FLOW 

Consider  a  uniform,  spacewise  and  timewise,  supersonic/ 
hypersonic  flow  of  an  inviscid,  perfect  gas,  with  constant 
specific  heats,  past  a  formation  of  two  wedges  at  design 
condition  (Figure  3.1a).   The  wedges  are  assumed  to  be 
oscillating  in  pitch,  with  small  amplitude  and  frequency, 
in  an  independent  fashion  and  the  bow  shocks  are  assumed  to 
be  attached  to  both  bodies  at  all  times.   The  flow  quantities 
over  the  upper  and  lower  sides  of  the  second  two-dimensional 
wedge  (in  Regions  C  and  D)  are  required.   It  is  clear  that, 
depending  on  the  difference  of  the  semi-vertex  angles  of  the 
two  wedges ,  either 

a)  both  sides  of  the  second  wedge  are  compression  sides 
(Figure  3.1a)  or, 

b)  one  side  of  the  second  wedge  is  a  compression  side  and 

2 
the  other  side  is  an  expansion  side  (Figure  3.1b) . 

These  compression  and  expansion  side  problems  are  con- 
sidered in  subsequent  subsections. 

In  subsections  A  and  B  the  flow  quantities  over  the  com- 
pression and  expansion  sides  of  an  oscillating  wedge  will  be 
given.   The  upstream  flow  (in  Region  B) ,  is  assumed  to  be 


This  problem,  suggested  by  Prof.  Platzer,  is  of  interest 
in  high-speed  turbomachinery  aerodynamics-. 
2 
The.  possibility-  of  the  second  side  being  neither  a 
compression  nor  an-  expansion  side  is  not  excluded.   In  this 
case  though  region  C  becomes  an  extension  of  region  B,  in 
which  region   the  flow  quantities  are  assumed  to  be  known. 
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shock 


edge    H 


shock 


Fig.   3.1a       Formation       of     two     wedges.     Both      sides,    of 
wedge    XE       are      compression       sides. 
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Fig.    3. lb       Fonmafion     of      two     wedges.      We-dge     II     has     one 
expansion    and     one    compression      side. 
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oscillatory  with  circular  frequency  u>  and  flow  quantities 
of  the  form 


where : 


u  =   uQ(l  +  eU)  (3.1a) 


v  =   uq  £  V  (3.1b) 


p   =  p  (1  +  eyMP)  (3.1c) 

o        o 

p   =   P0(l  +  £MQR)  (3. Id) 


u,  v  are  the  x-  and  y-velocity  components  relative 
to  the  coordinate  system  shown  (Figure  3.1a); 


A  ikt 
e  =  ee     with  e :   a  small  parameter  characteristic  of 

the  deviation  of  the  upstream  flow  quantities 

from  their  average  value; 

k:  the  so-called  reduced  frequency  defined  by 

k  =  a£  ; 

u 
o 

I:    a  characteristic  length  of  the  flow; 

u  ,p  ,  p   and  M  are  the  mean  or  average  velocity,  pressure, 
density  and  Mach  number  of  the  upstream  flow; 

U,V,P  and  R:  are  time  independent  quantities  considered  to 
be  known  functions  of  the  non-dimensional 
spatial  coordinates  x  and  y. 


We  will  restrictthe  analysis  to  small  reduced  frequencies 
k  and  will  assume  that  these  quantities  have  the  form 

0  =  U(0)  +  (ik)U(1)  =  u<0)+(ik)  (u^'x+uf 'y+u*1')        (3.2a) 

V  =  V(0)  +  (ik)V(1)  =  v<0)+(ik)  [vf'xtvf1^1')        (3.2b) 
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p  =  p'0'  +  (ik)P(1'  =  p<°»  +  (ik,  Ipf'^jVp]1'!        (3.2  =  ) 
R  =  R(0»  +  (ik)R(1)  =  r<0)  +  (ik)  Cr'^x+rf'y+r'1')       (3. 2d) 

where  the  coefficients  u^  '  ,  u.;   ,  ...  are  known  constants. 

It  should  be  noted  that  the  above  forms  contain  as 
special  cases  the  closed  form  solutions  found  in  Sections 
III.A.l  and  III. A. 2  for  the  flow  over  the  upper  and  lower 
sides  of  an  oscillating  flat  plate  at  an  arbitrary  angle  of 
attack.   Thus  the  results  given  in  the  following  subsections 
A  and  B  hold  true  for  upstream  oscillatory  fields  of  the 
general  form  (3.2)  and  are  not  limited  to  oscillatory  fields 
created  by  oscillating  flat  plates/wedges.   Similarly  the 
solution  given  in  subsection  A  contains  as  a  special  case 
the  solution  to  the  problem  of  a  stationary  wedge  in  an 
oscillatory,  uniform  spacewise,  hypersonic  free  stream,  which 
was  studied  in  [Ref.  9].   Furthermore  the  solution  presented 
here  is  exact  and  holds  for  the  complete  supersonic/hypersonic 
speed  range. 

In  subsection  C  an  alternate  approach  to  the  expansion 
side  problem  is  suggested.   In  this  approach  the  upstream 
flow  is  assumed  to  be  oscillatory  and,  for  simplicity,  uni- 
form spacewise .   Two  boundary  value  problems  with  linear 
equations  and  boundary  conditions  are  formulated  and  a  closed 
form  solution  of  the  problem  is  sought. 
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A.   UNSTEADY,  NON-UNIFORM  FLOW  PAST  AN  OSCILLATING  TWO- 
DIMENSIONAL  WEDGE — COMPRESSION  SIDE 

Consider  a  wedge  (Body  I)  ,  at  design  condition,  osci- 
llating with    small  amplitude  and  frequency  in  supersonic/ 
hypersonic  flow  (Figure  3.2) .   This  flow  problem  was  studied 
in  Section  III.A.l  and  the  flow  field  quantities  u_ ,  v. ,  p _ , 

B    B    B 

p  ,  ...  in  region  B  (expressed  relative  to  axes  O'x'y' )  were 
completely  determined.   Assume  now  that  a  second  wedge 
(Body  II) ,  located  entirely  in  Region  B,  is  oscillating  also 
with  small  amplitude  and  frequency.   The  flow  field  quanti- 
ties u  ,  v  ,  p  ,  p  ,  ...  in  Region  C  are  required. 

V—  (w  \^  C_* 

To  solve  the  problem  we  adopt  the  following  procedure. 

a)  Express  all  flow  field  quantities  relative  to  the 
coordinate  axes  Oxy  attached  to  body  II  in  its  mean  position. 

b)  Assume  that  both  bodies  I  and  II  are  stationary  and 
find  the  reference  steady  flow  quantities  in  Region  C. 

c)  Assume  that  body  I  is  oscillating  while  body  II  is 
kept  stationary  and  superimpose  to  the  steady  flow  quanti- 
ties in  region  C  perturbation  quantities  due  to  the  oscilla- 
tion of  body  I . 

d)  Assume  next  that  body  II  also  is  oscillating  and 
superimpose  to  the  aleady  perturbed  flow  quantities  in  Region 
C  new  perturbation  quantities  due  to  the  oscillation  of  body 
II. 


Flow  quantities  in  region  D  can  be  found  in  a  completely 
analogous  way  (by  simply  letting  9,  =  9,  +  92  in  the  solution) 
since  the  flows  are  assumed  independent. 
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Tig. 3. 2  Illustration       of       the     flow     quantities. 
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The  procedure  described  above  effectively  breaks  down 
the  flow  problem  in  hand  (namely,  the  calculation  of  the 
flow  quantities  in  Region  C  when  both  wedges  are  oscillating) 
into  a  sequence  of  three  separate  problems  which  may  be  solved 
successively  to  give  the  final  results . 

The  solution  method,  which  is  based  on  Professor  Hui ' s 

theory,  is  presented  in  Sections  IV. A. 2  and  IV. A. 3.   The 

final  results  are  given  in  Section  IV. A. 4.   In  the  following 

section  IV.A.l,  the  problem  is  formulated  and  the  flow  field 

quantities  in  Region  B  are  expressed  relative  to  axes  Oxy 

attached  to  body  II  in  its  mean  position. 

1.   Formulation  of  Problem — Expressing  Upstream  Quantities 
Relative  to  Coordinate  System  Attached  to  Second  Body 

Consider  steady  uniform  supersonic/hypersonic  flow 

past  the  formation  of  the  wedges  described  above  (Figure  3.3). 

Assume  that  the  wedges  have  chord  lengths  £,  and  T~    and  are 

oscillating  with  circular  frequencies  uj,  and  a^.   Cartesian 

coordinate  systems  O'x'y'  and  Oxy  are  attached  to  the  bodies 

with  origins  placed  at  the  wedge  apexes  and  axes  Q'x1,  Ox 

along  the  mean  positions  of  the  upper  surfaces.   Let  the 

steady  flow  quantities  in  Regions  B  and  C  (stationary  wedges) 

be  given  by  u   ,  p   ,  p   and  u  ,  p  ,  p  ,  respectively.   Also 
3       2      o_ '   o  '  ^o_,      0    0*0 

ODD 

denote  by  t  the  time  variable  and  by  u  ,  v  ,  p   and  p   the 
perturbed  quantities  in  Region  B  expressed  relative  to  axes 
O'x'y' .   The  departure  of  the  perturbed  flow  from  the  steady 
flow  in  Region  B  is  characterized  by  the  small  quantity  z~ 
introduced  below  and  the  reduced  frequency  associated  with 
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U-.    M£ 


Fig.    3.3  Illustration     of      coordinate       systems. 
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the  oscillation  of  the  first  wedge  is  defined,  as  usual,  by 

k,  =  — .   Similarly  the  reduced  frequency  associated  with 

°B                                            (ja25'2 
the  oscillation  of  the  second  wedge  is  defined  by  k~  =  

o 

Non-dimensional  time  and  lengths  are  introduced  associated 

with  bodies  I  and  II  as  follows, 


t]_   =   t/£1#    t2   =  t/l2,        x'   =   x*A1,   x  =  x/i2 


y'  =  y'A^  y  =  yA2 

We  assume  that  the  perturbed  quantities  in  Region  B 
have  the  form 

u_   =  u   (1  +  e.U)  (3.3a) 

B      oB      1 

vn   =   e,  u   V  (3.3b) 

B       1   oB 

-   =   p   (1  +  e,YM   P)  (3.3c) 

PB       °B       X   °B 

p_   =   p   (1  +  e.M^  R)  (3.3d) 

B      °B       L   °B 

ik.t 
where  e.  =  e,e     and 

U   =  u^0)  +(ikx)  (u^1}x'  +u^1}y'  +u^1})  (3.4a) 

V  =  v^0)  +(ikx)  (v^x'  +v^1)y'  +v^1))  (3.4b) 

P   =  p<0)  +(ik1)(p{1)x'+P^1)y+pf))  (3.4c) 

-  r<0)  +(ik1)(ri1)x'+r«1)y+r«1))  <3.4d) 


R 
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with  coefficients  u3   ,  uj   ,  ...  assumed  to  be  given.   We 
recall  that  in  Sections  III.A.l  and  III. A. 2  the  problems 
of  uniform,  steady  supersonic/hypersonic  flow,  past  harmonically 
oscillating  wedges/flat  plates  were  studied  and  the  perturbed 
quantities  behind  the  shocks/expansion  fans  were,  in  both 
cases,  found  to  be  of  the  form  assumed  above  by  relations  (3.3) 
and  (3.4).   Thus  the  analysis  that  follows  does  not  distinguish 
whether  the  disturbed  flow  field  in  Region  B  has  originated 
from  an  oscillating  wedge  or  flat  plate  over  its  expansion 
and  compression  sides  and  the  appropriate  coefficients  and 
parameters  should  be  used  to  make  the  distinction  for  the 
specific  case  considered. 

Next  we  express  the  quantities  given  by  (3.3)  and 
(3.4)  relative  to  the  coordinate  system  0  xy  attached  to 
body  II.   Parallel  transformation  and  rotation  of  axes  gives 
for  the  general  quantity 


<°>  +  (^^'^.fv*.;11) 


M  =  m_ 


where  symbols  M  and  m  stand  for  capital  and  lower  case  sym- 
bols U,  V,  ...  and  u,  v,  ...  respectively. 


M  =  m<0)  +  (ik.)C'  x  +  C '  y  + C '  )  (3.5) 

3        1  nu      m^    m~ 


where 


a   =  a  cos  91  +  31  sin  9^   3Q  =  3X  cos  Q1   -  c^  sin  61, 
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9d     ■      92    "   91 


\     "      (m^1'  cos9d+m^1)  sined)^2/£1 


Cm2     m      <m2l)  =°sed-m^1)   sined)£2/F1 

c;3  ■  n,31)  +  (aoml1)+Boml1)»vi1 


Combining  relations  (3.3)  and  (3.5)  we  express  the  perturbed 
quantities  pB  and  p"  relative  to  axes  Oxy  in  the  following 


form 


p   =  p   (1+e  [C   +  (ik  )(C  x+C   y+C   )]}            (3.6a) 

B  °B      1      PQ     1      Pl    P2    P3 

p.   =  p   (1  +e,  [C   +(ilO  (C   x+C   y+C   )]}           (3.6b) 

B  °B     1      ro    1   rl   r2   r3 

with 


C    =  y   M   p*0),     C  =  M   r'0) 
Po        °B   3        ro      °B   3 

C    =  v  M  C  ,      C  =  M   C    (j  =  1,2,3) 
Pj         °B  Pj         r j       °B   r j 

To  find  the  perturbed  velocity  components  u_ ,  v_  we  take 

their  components  along  the  new  axes  and  add  them  to  get 


u_   =  u   [cos  9  ,  +  £,  (U  cos  9  ,  +  V  sin  9  ,)  ] 
B      o  d   1       d        a 

a 


vD      =      u      [-  sin  6,  +en  (-U  sin  9  ,  +V  cos  9  ,)  ] 
B  o„  d        1  d  d 


118 


and  finally  using  relations  (3.5)  we  have 


u    =   u   (cos  ed  +  6l[C   +(ik  )(C   x+C   y+C   )]}   (3.6c) 
B  o         1     2     3 

^B      =      Uo    {~sinQd  +el[Cv  +(iki)(C      x+C      y+C      )]}       (3.6d) 
B  vo  1         v2  3 

where : 


C    =  u<0) cose, +v<0) sine,,   C    =  -u<0) sine  ,+v^0) cose^ 
uQ       3      d3      dv        3      d3      d 

Cu.   =   Cu.COsed+Cv.Sined   (j  =  1'2'2) 

Cv   =   -C^  sined+C^  cose   (j  =  1,2,3) 
3         j         J 

Relations  (3.6)  give  the  complete  set  of  flow  quantities  in 

Region  B  expressed  relative  to  coordinate  system  0  xy. 

2.   Perturbed  Equations  of  Motion 

In  this  section  we  derive  the  perturbed  equations  of 

motion  in  Region  C,  when  both  wedges  are  oscillating  (Figure 

3.2).   To  find  the  form  of  the  perturbed  quantities  in  this 

region,  relative  to  axes  0  xy ,  we  proceed  as  follows. 

First  we  assume  that  both  bodies  are  stationary  and 

denote  the  (reference)  steady  flow  quantities  in  Region  C 

by  u  ,  p  ,  p  ,  M   (Figure  3.4a) .   Obviously  v  =  0  for  the 
-1      o   o   o   o  o 

coordinate  system  chosen. 

Next  we  assume  that  body  I  is  oscillating  while 
body  II  is  kept  stationary.   Due  to  the  oscillation  of  body 
I  we  superimpose  to  the  field  quantities  perturbations  of 
the  form  (Figure  3.4b), 

119 


G.=  -  y+x tonvf  =0 


Body    II 


Fig.  3. 4a  Flow      quantities     in       Region      C.        B  o  t.h 

wedges      are      stationary.  - 


•  u 
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g.   3.4b  Flow      quantities    in      Region    C.       Body 

oscillating       w"hile    Bo-dy    II     is     stationary. 
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ik   t 
£1U1     =      £le      1   luoUl     =      £luoUl  (3*7a) 

ik   t 
=lvl     =     £le     1   luoVl     =     eluoVl  (3-7b) 

ik  t 
*1*1     =     "le      X   Wl     =     £lpo™oPl  (3-7c) 

e-pn       =      e.e      X    1p   M   R.       =      e.p   M   R,  (3.7d) 

11  1  o  o   1  1  o   o   1 


Finally  we  assume  that  both  bodies  are  oscillating. 
Due  to  the  oscillation  of  body  II  we  superimpose  to  the  field 
quantities  perturbations  of  the  form  (Figure  3.4c) 

ik  t 

e,u.   =   e0e   2  2u  U.   =   e-u  U,  (3.8a) 

21       2        o2       2o2 

e0vn   =   e0e  *      u  V_   =   e0u  V0  (3.8b) 

21       2        o2    .   2  o  2 


ik2t2 


£2P1   =   £2e      PoYMoP2   =   £2PoYMoP2 

ik  t 
£2P1   =   £2e   2  2poMoR2   =   £2poMoR2 


(3.8c) 


ik  t 

4  ^0  M  R„   =   e„p  M  R„  (3.8d) 


The  flow  field  quantities  in  Region  C  for  the  last  case 

2  2 

become    (terms  of  order  e-^    '    ei    '    £2     are   ne<?lected) 


u      =      u    (1   +   e.U.    +   £0U0)  (3.9a) 

O  11  2     2 


v      =      u    (e.V.    +   £0V„)  (3.9b) 

Oil  22 


p      =      p    [1  +  YMo(£lP1+e2P2)]  (3.9c) 
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♦    fc  (  *-h  cos62  )=  0 


Fig.    3.4c         Flow       quantities     i 
bodies      oscillating. 


n       Re  gion     c   .       Both 
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p   =   pQ[l  +  Mq(£1R1  +e2R2)]  (3.9d) 


where  the  quantities  U,,  U  ,  v.,  ...  are  functions  of  the 
non-dimensional  variables  x  and  y.   For  small  k   and  k 
we  may  expand  these  quantities  as  power  series  in  (ik)  of 
the  general  form 

X.   =   xf0)  +  (ik.)XU)  +  ...   (j  =  1,2)      (3.10) 
3  3  3      3  J 

with  X  denoting  U,  V,  P  and  R. 

To  simplify  we  assume  that  the  characteristic  lengths 
J,  and  £_  are  equal  and  we  put  (3.9)  into  the  governing 
equations  of  motion  (2.1)  .   Using  (3.10)  and  equating  the 
terms  of  the  same  order  of  (ik)  in  each  of  the  resulting 
equations  we  obtain 


el(V<°>  +V<°>  +MoR<°>,  +  e2<U<°>  +  V<°>  +MoR^')  =  0  O.lld, 

^(ik.Xu  u{0'+uoU^»)+e2(ik2)uo(U^)+U^») 
B 

=   -Vel(ikl)Pl°>+£2(ik2)P2x)!/Mo         (3-12a) 
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£l(ikl)(Vl0,+U°V^,)+e2(ik2)uo(V2°>+V2x)) 

0 

=   -uo[e:1(ik1)p{°)+e2(ik2)P^)]/Mo       (3.12b) 
1   1   o_.   1    1     o   lx   lx 


B 


+  e,(ik9)u  (P^0)-R^0)+P#[1)-R^))  =  0      (3.12c) 
2    2   O   2     2     2x    2x 


e.Uk-JCu  MR.(0)+u  (U.(1)+V,(1)+M  R.(1))  ] 
1    1    o      o  1     o   lx    ly    o  lx 


iy 

+  £2(ik2)uo(MoRr)+U2x)+Viy)+MoR2x))  =  °  (3'12d) 

Equations  (3.11  and  (3.12)  are  the  zeroth  and  first-order 
equations  for  our  problem.   For  the  case  of  a  single  oscillating 
body  they  reduce  to  equations  (2.9,  2.10)  given  in  Section 
III.A.l,  by  letting  e   =  0 ,  U   =  V   =  P   =  R  =  0  (body  I 
considered  missing)  and  dropping  subscripts  B  and2  in  re- 
sulting expressions. 

To  solve  for  the  sixteen  unknowns  contained  in 
Systems  (3.11)  and  (3.12)  we  may  proceed  in  two  steps.   We 
may  first  find  the  eight  unknowns  associated  with  the  oscilla- 
tion of  body  I.    To  do  this  we  will  assume  that  body  II  is 
not  oscillating  (e_  =  0)  and  solve  successively  the  above 
systems  of  equations  subject  to  appropriately  formulated 

boundary  conditions  for  this  case.   Next,  we  will  find  the 

2 
eight  unknowns  associated  with  the  oscillation  of  body  II. 


Unknowns  with  subscript  1 . 

2 

Unknowns  with  subscript  2. 
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To  do  this  we  will  assume  that  body  II  is  oscillating  while 
body  I  is  present  but  stationary  (e,  =  0)  and  solve  successively 
the  resulting  systems  of  equations  subject  to  boundary  con- 
ditions appropriate  for  this  case.   The  second  problem  has 
been  treated  in  Section  II. A  and  we  will  not  repeat  the  solu- 
tion for  this  case.   We  note  though  that  the  problem  should 
be  appropriately  reformulated  to  account  for  the  different 
direction  and  magnitude  of  the  approaching  steady  flow 
velocity  as  indicated  by  the  ratio  of  steady  velocities 

u  /u   appearing  in  equations  (3.12). 
s 

3 .   Boundary  Conditions 

For  both  wedges  oscillating,  let  the  equation  of  the 
surface  of  body  II  be  given  by  (Figures  3.4b,  3.4c), 

S   =   -y  +  e2(x  -  h2cos92)   =   0 

and   the   equation   of    the   shock   attached    to   body   II    be   given 
by    (Figures  3,4), 

G      =      -y   +   x  tan  $   +   £-,Q,  (x)    +   e2^ix^      =      ° 

where  Q. (x)  and  Q,, (x)  are  unknown  quantities  associated  with 
the  oscillations  of  bodies  I  and  II  respectively  and  are  to 
be  found  as  part  of  the  solution. 

We  will  formulate  the  boundary  conditions  for  the 
case  of  stationary  wedge  II  (e_  =  0) . 

The  flow  tangency  condition  (2.2)  with  V  =  u0^1+£iui' 
e.V..}  gives 
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v{°}   =0      At   y   =   0  (3.13a) 

V^1)   =0      At  y   =   0  (3.13b) 


To  find  the  boundary  conditions  across  the  shock  we 
use  the  Rankine-Hugoniot  conditions.   The  procedure  is  given 
in  Appendix  D  and  the  resulting  zeroth  and  first-order  condi- 
tions, after  linearization  and  use  of  equation  (3.10),  have 
the  following  general  form 

Y1°}   =  KY  Ql(0)  +  KY  At  y  =  xtan$   (3.14a) 

x 

Y<1}   =   K    x+K    y  +K  +K   Q'(1)  +K   Q<0) 
xx     yy         x         o 

At  y  =x  tan  (J)    (3.14b) 

where  Y  stands  for  U,  V,  P  and  R  and  the  coefficients  K  , 
K   ,  K    /  . . .  are  known  constants  (functions  of  geometry 

X      XX 

and  steady  flow  quantities)  given  in  Appendix  D. 
4 .   Solution 

Two  boundary  value  problems  have  been  set  up  and 
will  be  solved  successively.   The  zeroth-order  equations  1,3 .11) 
with  £p  =  0  and  boundary  conditions  (3.13a)  and  (3.14a)  con- 
stitute the  first  problem,  which  will  give  the  steady  flow 
quantities  behind  the  shock,  for  the  case  of  stationary 
wedges  in  uniform,  steady  supersonic/hypersonic  flow.   The 
first-order  equations  (3.12)  with  e2  =  0  and  boundary  condi- 
tions (3.13b)  and  (3.14b)  constitute  the  second  problem.   Since 
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both  problems  are  linear  we  assume  that  their  solutions 

are  linear  combinations  of  the  non-dimensional  coordinates  x 

and  y. 

For    the   first  problem  we   assume   a   solution  of   the 
form 

U<°>      =      u<°»x  +   u<°'y   ♦  u<°> 


V<°>      =      v<0)x  +  v<0)y   +  v<0) 


R{°>   =  r<°»x  +  r'O'y  +  ,») 


1       Pl      p2   y   P3 


n(0)       (0)      (0) 
Ql       ql      q2 


Substitution  of  these  values  in  the  zeroth  order  equations 
and  boundary  conditions  gives  the  following  solution. 


v{0)   =   0  (3.15b) 

>r  -  pjo)  =  kpx^o) 

Q(0)   =    (0)x  .  _JV_X  (3.l5e) 

11         "Vx 
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For  the  second  problem  we  assume  a  solution  of  the 
form 


u'1'   =  u»>»  +  uf  >y  +  u'1' 


v'1'   =  v^'x  +  vf'y  +  v'1' 


R*1'   =  r^'x  +  r<X>y  +  r*1' 


p'1'   -  p<»*  +  p«»y  +  p<» 


.11)      „(D,2    „U)V    _(D 

Q1     =   q1   x   +  q2   x  +  q3 


Substitution  in  the  first-order  equations  and  boundary  con- 
ditions gives  the  solution 


U*1'  =  u^'x  ♦  ««»,  +  u'1'  (3.16a) 

V^11  =  v^'y  (3.16b) 

R*1'  =  r^'x  +  r^'y  +  r'1'  (3.16o) 

P*1'  =  pjW,  +P*1'  (3.16d) 


with 


V2X)   =   (KVxql  +  WV)/tan  * 


P-.(1)   =   K„  q,  +  W_ 
1         Px^l     P 


129 


P(1)       -      K      a(1)     +    K 
P3  ~      KPxq2         +    KP 


■J"      -      -^'/M^uf'u      /uo, 

B 

U21}       =       (~U11)     +   KUxql    +   V/tan    * 


u(1)       =      K      a(1)    +   K. 
U3  KUxq2        +   *U 


r«"      -     Pi1'    +u  o(Pi0»-r<°»,/uo 

0 

r^1)   =   (-r:[1)  +  KRxq1   +   WR) /tan  <j> 


r(1)   -   K   a(1)+  k 
r3     "   KRxq2   +  KR 


q21}   -   -V*7x 

a    -   2a(1>       f^n^   n<°>M  \         (MQ~1)WP 

ql   "   2qi     =   [U-(U3    ~P3   Mo) M 

o  o 

W      (M2-1)K„     K.y 
v  ]/[   S  ,   Px  +  /x  J 


tan  <J)       M       tan  4 
o 


(0) 

Y       Y 

xx     yy  o 


W„   =   K„    +  Ky    tan  $  +  K^  q  * 


with  Y  standing  for  U,  V,  P  and  R. 

B.   UNSTEADY  NON-UNIFORM  FLOW  PAST  AN  OSCILLATING  TWO- 
DIMENSIONAL  WEDGE — EXPANSION  SIDE 

Consider  a  wedge  (Body  I),  at  design  condition,  oscillating 

with  small  amplitude  and  frequency,  in  supersonic/hypersonic , 

steady,  uniform  flow  (Figure  3.5).   This  flow  problem  was 

studied  in  Section  III.A.l  and  the  flow  field  quantities  in 
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Fig. 3.5  I   llustration      of      the     flow     quantities      for     the 

two        wedge        problem      .     expansion     case       shown. 
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Region  B,  expressed  relative  to  axes  O'x'y'  were  completely 
determined.   Assume  now  that  a  second  wedge  (Body  II) , 
loaated  entirely  in  Region  B,  is  also  oscillating  with 
small  amplitude  and  frequency.   Let  the  semi -vertex  angles 
of  the  wedges  be  If.  and  V2   respectively  and  assume  that  their 
difference  V2   -¥-  is  negative.   In  this  case  the  upper  side 
of  the  second  wedge  (Region  D)  becomes  an  expansion  side  and 
the  exact  solution  given  in  the  previous  subsection  does  not 
app  ly . 

To  solve  the  problem  we  adopt,  as  before,  the  following 
procedure. 

a)  We  assume  that  both  bodies  are  stationary  and  find 
the  reference  flow  quantities  inside  the  expansion  fan 
(Region  C)  and  over  the  upper  surface  of  the  wedge  (Region 
D)  . 

b)  We  assume  that  body  I  is  oscillating,  while  body  II 
is  kept  stationary  and  superimpose  to  the  steady  flow  quan- 
tities in  Regions  C  and  D  perturbation  quantities  due  to  the 
oscillation  of  body  I. 

c)  We  finally  assume  that  body  II  is  also  oscillating  and 
superimpose  to  the  already  perturbed  flow  quantities  in 
Region  C  new  perturbation  quantities  due  to  the  oscillation 
of  body  II.    The  problem  is  thus  solved  in  three  successive 


Disturbances  due  to  the  oscillation  of  body  II  propagate 
along  the  characteristics  and  are  not  expected  to  influence 
Region  C  except  in  the  neighborhood  of  line  OC . 


13  2 


steps.   In  what  follows  we  will  restrict  the  analysis  to 

the  case  of  a  stationary  wedge  in  oscillatory  flow.   The 

solution  can  be  extended  to  the  case  of  an  oscillating  wedge 

by  following  the  procedure  described  in  Section  III. A. 2. 

In  subsections  1  and  2  the  problem  is  formulated  and  the 

method  of  solution  is  presented.   This  method  calculates  the 

flow  quantities  in  the  expansion  fan,  in  a  sweeping  fashion 

from  Region  B  to  Region  D,  along   rays  9  =  constant,  with 

8  increasing  in  small  steps  from  0 (M   )  to  9 (M)  (Figure  3.6). 

°B 
The  procedure  involves  the  repeated  application  of  two  basic 

steps  and  is  described  in  subsection  3.   The  solution  is 

approximate  and  becomes  more  exact  in  the  limit  as  the  number 

of  iterations  increases. 

1.   Formulation  of  Problem — Expressing  Upstream  Quantities 
Relative  to  Coordinate  System  Attached  to  Second  Body 

Consider  steady,  uniform,  supersonic/hypersonic  flow 

past  the  formation  of  the  wedges  described  above  (Figures, 

3.5,  3.6).   Assume  that  the  wedges  have  chord  lengths  &,  and 

I"2  and  are  oscillating  with  circular  frequencies  go,  and  cu2. 

Cartesian  coordinate  systems  O'x'y'  and  Oxy  are  attached  to 

the  bodies  with  origins  placed  at  the  wedge  apexes  and  axes 

O'x^Ox  along  the  mean  positions  of  the  upper  surfaces.   Let 

the  steady  flow  quantities  in  Region  B  be  given,  relative 

to  axes  O'x'y1 ,  by  u   ,  p   ,  p   ,  M    ...  and  in  Region  D, 

_°B    °B    °B    °B 
relative  to  axes  Oxy,  byu,  p,  p,  M,  ...  .   Let  also  the 

steady  quantities  along  a  ray  9  =  constant  in  the  Prandtl- 

Meyer  expansion  fan  be  given,  relative  to  axes  Oxy,  by 
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U*B    ■  M*B 


Vr 


shock 


Fig. 3. 6         Polar     coordinate        sysfem       attached       to 
body     I  I   . 
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ue'  ve '  pe'  pe'  Me  •••  where 

ue  =  vr  cos  (r"e)  +  vesin(r-e) 
vQ  =  vr  sin  (r-e)  -v  cos  (r-e) 

where  the  angle  r  and  the  r-  and  8-  velocity  components 

V  and  VQ  were  defined  in  Section  II.A.4.b. 
r      o 

We  assume,  as  in  Section  IV. A. 1,  that  the  perturbed 
quantities  in  Region  B  have  the  form  given  by  equations  (3.3) 
and  (3.4)  and  note  as  before  that  this  form  contains  as 
special  cases  the  form  of  the  perturbed  quantities  behind 
oscillating  flat  plates/wedges .   To  express  the  quantities 
given  by  (3.3)  and  (3.4)  relative  to  the  coordinate  system 
0  x  y  attached  to  body  II  we  use  relations  (3.6)  with  9,  =  ^"^l 
2.   Method  of  Solution 

Assume  that  the  change  in  flow  direction  over  the 
corner  (Figure  3.7a)  is  obtained  in  n  steps  (Figure  3.7b) . 
Assume  that  the  step  changes  in  flow  direction  are  all  equal, 
i.e. , 


a,  =  an   =    ...  =  a  =  — ,     a  -  ¥,-¥., 
12  n   n  1  2. 


and  also  note  that  each  of  the  line  segments  DkDk+2_  (k  ■  1/2, 
...,n-l)  extends  to  infinity.   In  the  second  case  the  single 
expansion  fan  has  been  replaced  by  n  smaller  expansion  fans 
and  the  two  problems  are,  physically,  completely  equivalent. 
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REGION     B 


Fig.  3.7         Change      i 


n     flow     direction.  (a)    Over      a 


single    corner         (b )       Over      a       series       of        smaller 


c  orners 


136 


Next  we  assume  that  the  change  in  flow  direction  over 
the  step  corners  is  achieved  through  expansion  discontinui- 
ties, similar  in  nature  to  hypothetical  expansion  shocks, 
which  we  will  call  expansion  fronts  or,  simply,  fronts. 
Suitably  spaced  lines  0B1 ,  0B2,  . . . ,  OB   in  (Figure  3.7a) 
correspond  to  the  mean  positions  of  the  fronts  F,,  F,,  ..., 
F   emanating  from  corners  D,  ,  D2,  ...,  D   respectively. 
These  lines  divide  Region  C  in  subregions  B  ,  B,  ,  B2,  ..., 
Bn_3/  Bn  E  D*   The  first  subregion  B   is  separated  from 
Region  B  by  line  OB  =  OB  and  from  Region  B,  by  front  F, . 
Subregions  B,   (k  =  l,2,...,n-l)  correspond  to  line  segments 
D,  D,  -  and  are  separated  from  their  adjacent  subregions  from 
the  left  by  fronts  F,  and  from  the  right  by  fronts  F,  , .   The 
last  subregion  B   ED  corresponds  to  the  line  segment 
D  D   -  E  D  D  (actual  upper  surface  of  wedge)  and  is  separated 
from  Region  C  by  front  F  .   The  flow  quantities  in  this  last 
subregion  are  required.   We  intend  to  solve  for  the  unknown 
quantities  in  this  region  by  first  obtaining  the  unknown 
quantities  in  Regions  B, ,  B-,  . . . ,  B   ,  successively. 

We  observe  that  the  original  flow  problem  has  been 
replaced  by  a  set  of  n  identical  flow  problems  each  of  which 
involves  the  same  step  change  in  the  direction  of  flow.   We 
may  further  observe  that  the  assumption  of  "jump"  step  changes 
in  the  direction  of  flow,  via  the  so-called  fronts,  is  simi- 
lar to  the  "jump"  step  changes  in  the  direction  of  flow  via 
ordinary  compression  shocks.   Thus  the  two  problems  shown  in 
(Figures  3.8a  and  3.8b)  are  mathematically  identical  if  the 
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Figure  3.8.   Duality  of  problems 
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common  term  discontinuity  is  used  to  denote  "jump"  changes 
in  flow  direction  and  quantities  across  it. 

In  Section  IV. A  the  solution  to  the  flow  problem 
shown  in  (Figure  3.8a)  was  obtained.   It  was  assumed  there 
that  the  approaching  flow  was  oscillatory  of  the  form  given 
by  equations  (3.3)  and  (3.4)  and  the  flow  quantities  in  Region 
B,  were  found  relative  to  the  coordinate  system  D  x,  y..   We 
will  see  that  the  approaching  flow  for  the  dual  problem  shown 
in  (Figure  3.8b)  is  of  the  same  form  and  thus  the  solution 
given  in  Section  IV. A. 4  applies  to  the  dual  problem  also, 
provided  that  the  appropriate  geometrical  data  and  steady 
flow  quantities  are  introduced  for  this  case.   The  quantities 
found  in  Region  B,  will  be  expressed  relative  to  the  next 
coordinate  system  D.  ,  Xv+tYv+i  an<^  will  be  used  as  input 
quantities  to  an  identical  problem  to  give  the  solution  for 
the  flow  quantities  in  the  next  Region  B,  ,  .   The  procedure 
can  be  repeated  till  the  required  quantities  in  the  last 

region  D  =   D  are  found. 

r     n 

Physically  the  assumption  of  a  "jump"  expansion  dis- 
continuity is  not  accepted  since  such  a  finite  discontinuity 
would  lead,  as  explained  in  Section  II.A.3.b  to  a  decrease 
in  entropy.   The  procedure  described  above,  however,  should 
give  the  exact  solution  in  the  limit,  as  the  number  of  step 
changes  in  angle  is  increased. 
3 .   Solution  Procedure 

Consider  an  arbitrary  Region  B^  (k  =  1,2,..., n)  behind 
front  F.   (Figures  3,8b,  3.9).   Assume  that  both  bodies  are 
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£>  (  FAN    ANGLE  ) 


Figure  3.9.   Illustration  of  geometrical  data 
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stationary  and  let  the  reference  steady  flow  quantities  in 
this  region  be  u,  ,  v,  ,  p,  ,  p,  ,  M.  ,  ...  .   These  quantities 
can  be  found  from  the  Prandtl-Meyer  relations  given  in 
Section  II.A.4.b.   Expressed  relative  to  coordinate  system 
Ox,  y,  these  quantities  are 


P.  =  p^  [E(Mo  )/E(ML  )]Y/(Y-D                    (3.17a) 

K  Ob      Og  K 

p.  =  p„  [E(M   )/E(M,  )]1/(Y~1)                    (3.17b) 

k  oB     oB     Tc 

uk  =  Vr     COS  Ck   +  V9     sin5k                                                      (3.17c) 

K.  K 

vk  =  Vrksin^k  "  vekcos?k                                           (3*17d) 

Mk  =  [  (tan  X9k/X)2  +1]0,5                      (3.17e) 


where  p   ,  p   ,  M   ,  ...  are  the  reference  steady  flow 

°B    °B    °B 
quantities  in  Region  B  and 


X  =       (^-)°*5,     E(M)   =   l  +  (Y-l)M2/2 

9,       -      6     +k  i,  0         =      6(M      )       =      tan"1[X(M2    -1)°'5]/X 

k  o  n  o  oQ  oB 


V  =      c   sin    X9,,  VQ        =      X    c   cos    X9 

rk  k  9k 


c      .      0      [1  +  2         ,0.5 

°B  (Y-l)M^ 

°B 
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Assume  that  body  I  starts  oscillating  and  denote  the 

perturbed  quantities  in  Region  B,  expressed  relative  to 
coordinate  system  Ox  y  by 

u^.   =  uk(l  +  e^)  (3.18a) 

vk   =  ^(e^)  (3.18b) 

Pk   =   pk(1  +  elYMkPk)  (3.18c) 

Pk   =   pk(l  +  e^R^)  (3.18d) 


where  the  unknown  time-independent  quantities  U * ,  V '  ,  P/  , 
R/  may  be  expressed ,  for  small  amplitude  and  frequency  of 
oscillations,  as  power  series  in  (ik,)  of  the  form 

Yk  =  Yk(0)  +  <**].>  *£(1)  +  ••• 


with  Y'    standing   for  U',    V ,    P',    R1 . 

Assume    that,   with   body   I    oscillating,    the   equation 
of   the   front  F.     is   given  by 

*(Fk)      =      -^k   +  xk  tan*k  +elQk(xk}      =      °  (3.19) 


where 


k,  is  the  reduced  frequency  parameter  associated  with 
the  oscillations  of  body  I  and  should  not  be  confused  with 
the  integer  k . 
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^k   =   (T-e  -a)  +  (a-6)k/n,    r   =   a  +  ~  +  P (M   ) 


°B 


)   =   tan"1[X(M^  -1)°*5]/A  -  tan"1 (M2  -l)0*5 
B  °B  °B 


Assume  that  the  approaching  flow  from  Region  B, 
has  the  following  form,  with  flow  quantities  expressed  rela- 
tive to  system  Ox,y ,  , 

V-l   =   uk-l(cos  9d  +  £lUk-l}  <3-20a> 

*k-l   =  uk-l("sin  6d  +  £lVk-l}  (3'20b) 

Pk-1   =   pk-l(1  +  £lRk-l}  (3.20c) 

Pk-1   =   ?k-l(1  +  £lPk-l>  (3'20d> 

where  0   =  (m ^-^    ) /n   and  the  time  independent  quantities 

&v.i »  Vk-1'  ^k-l  and  Pk-1  are  assumed  to  be  known  quantities 
of  the  form 

Vl  "  CY  +  '"I'l^A+Vk^  J  (3-21) 

O  12        3 

with  Y  standing  for  U,  V,  R  and  P. 

We  observe  that  the  problem  formulated  above  for  the 
arbitrary  Region  B,  is  mathematically  equivalent  to  the 
problem  considered  in  Section  IV. A  since  the  approaching 
flow  has  the  same  form,  the  change  in  flow  quantities  occurs 
through  a  discontinuity  that  satisfies  the  same  equation  and 
therefore  the  same  Rankine-Hugoniot  boundary  conditions,  the 
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unknown  perturbed  flow  quantities  along  the  surface  have  the 
same  form  and  should  satisfy  the  same  flow  tangency  boundary 
condition  and  the  same  governing  equations  of  motion.   We 
may  therefore  borrow  from  Section  IV. A. 4  the  solution  for  our 
problem  as  given  by  equations  (3.15)  and  (3.16) .   As  a  final 
step  we  should  show  that  this  solution  expressed  relative  to 
the  coordinate  system  Ox.  ,,y.  ,,  (associated  with  the  problem 
that  is  to  be  considered  next  in  order  to  find  the  flow  quan- 
tities in  the  adjacent  Region  B.  ,)  is  of  the  form  assumed 
by  relations  (3.20)  and  (3.21).   We  see  that  this  is  indeed 
the  case  since  the  solution  given  by  equations  (3.22)  below 
is  expressed  relative  to  the  rotated  system  of  axes  OXfc+iYfc+i 
by  equations  (3.23). 

The  solution  given  by  equations  (3.9),  (3.15)  and 
(3.16)  is 

^   =   uk(l  +  £1[U]^(0)  +  (ik1)U]^(1)]}  (3.22a) 

-k  =  vi[vk(0)  +  (iv^1^  (3-22b) 

pk   =   pk(l  +  ^YM^P^0*  +  (ik^P^1*]}  (3.22c) 

Pk   =   Pk(l  +  e1Mk[B£<0)  +  (i*].)Rk(1)]}  (3*22d) 


with 


„.(0,  .  u(0)(   v,(o,  „  0f   ,.(0,  =  r(0)f 
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^(0)  -  pf1'  %a)  =  -i1'^*.^^". 


,(1)    _   „(1) 

k 


V'^'   =   V(1)v      R'(1)   -   r(1)x   +  r(1)v   +  r(1) 
V2   yk'    Rk      "   rl   Xk  +  r2   yk  +  r3   ' 


>im  ■  ^\'A" 


The  same  solution  expressed  relative  to  the  coordinate  axes 
0xk+lyk+l  is 

u.   =  u,  {cos9,+e,  [C   +(ik.  )(C   x,^,+C   yVj,+C„  ]}     (3.23a) 
k      k     d   1  u     1   u,  k+1  u^k+1  u0 

o         1        z        j 

v,       =      u.  {-sine  ,+en  [C      +(ik.)(C      x^.+C     y,  ^, +C      )]}       (3.23b) 
k  k  d     1     v  1       v,    u+1     v~Jk+l     v-, 

o  12  3 

Pk     -     Pk{1+el[S+(iV(S/*+l+S/k+l+V1}  <3-230) 

O  1  Z  J 

pk   =   pk{l  +.  [C  +(ik)(C  xk+1+C   yk+1+C   )]}       (3.23d) 

O  12  3 

with 


C  =      cos    9 ,u,      ,        C  =      cos    9 ,ui    ' , 

u  d   3      '  u-  d    3 

o  J 

C  =      cos    9-  (cos    9-ui    '    +   sin   9,ui        -   sin   6,vJ    }) 
u,  ddl  d2  d2 

C  =  cos  9,  (-sin  9,u-(1)  +  cos  9,u*    +  sin  9, v.;  J) 
u2  d       dl  d2  d2 

C  =   -sin  9Ju^0)  ,     C    =   -  sin  9  ,u^1) 
v  d  3        v.,  d  3 

o  -J 

C  =   sin  9d[-cos  9du[1)  -  sin  Q^"^  +  cos  9^  '] 
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C  =  sin   9,u1(1)    -   sin    6, cos    6,u^1)    +  cos29Jv^1) 
"-                         a   1                        a             a   2                         a    2 

C  -  Y^p'0*,              C           =      Y^p'1', 

cpj  =  TMk  cos  edP^  ,        cp2    =    YM,,  sin  edp»-> , 

C  =  Vcos   e^1'   -sin8dr^X>), 


cr2    ■    Vco"  vl11  +  sin  Vl"1' 


»d   ■    (VV/n 


We  summarize  the  procedure  that  should  be  followed 
to  find  the  flow  field  quantities  in  Region  B   =  D  below. 

(a)  Assume  that  the  flow  field  quantities  in  Region  B  have, 
relative  to  axes  O'x'y'  attached  to  body  I,  the  form  given 

by  relations  (3.3)  and  (3.4). 

(b)  Express  these  quantities  relative  to  axes  Ox,y1 ,  as 
described  in  Section  IV.A.l,  by  relations  of  the  form  given 
by  equations  (3.6),  with  9   =  ( ^-H^) /n. 

(c)  Use  these  quantities  as  upstream  quantities  for  the 
first  problem  considered  and  find  the  perturbed  quantities 
in  Region  B,,  as  described  in  Section  IV. A. 4.   The  solution 
has  the  form  given  by  equations  (3.2  2)  with  k  =  1  and  the 
reference  steady  flow  quantities  involved  are  given  by  equa- 
tions (3.17)  with  k  =  1. 
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(d)  Express  the  solution  relative  to  coordinate  system 
Ox2y2,  by  using  equations  (3.23)  with  8,  =  (^  -H'J/n. 

(e)  Repeat  step  (c)  to  find  perturbed  quantities  in  Region 
B2  and  repeat  step  (d)  to  express  them  relative  to  coordinate 
axes  Ox_y3. 

(f)  Repeat  step  (e)  till  the  flow  quantities  in  Region 

B   are  found.   These  quantities  are  the  required  quantities 
in  Region  D. 

The  procedure  given  is  well  suited  for  computer 
applications.   The  program  can  be  set  up  in  a  fairly  easy 
way  using  the  formulas  presented  and  the  number  of  iterations 
can  be  increased  to  the  accuracy  desired. 

C.   AN  ALTERNATE  APPROACH  TO  THE  EXPANSION  SIDE  PROBLEM 

In  the  last  section  the  flow  field  quantities  over  the 
expansion  side  of  the  wedge  were  found  by  a  series  of  itera- 
tive calculations.   The  same  technique  can  give  the  flow 
quantities  at  any  point  (r,8)  in  the  expansion  fan  region 
(Figure  3.10)  . 

In  this  section  another  approach  for  finding  the  flow 
quantities  in  the  expansion  fan  is  suggested.   The  approach 
consists,  as  before,  of  the  following  steps. 

a)  Expressing  the  governing  equations  of  motion  in  the 
expansion  fan  region  in  polar  coordinates. 

b)  Perturbing  the  equations  of  motion. 

c)  Forming  systems  of  equations  that  are  to  be  solved 
subject  to  appropriately  formulated  boundary  conditions. 
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Figure  3.10.   Illustration  of  polar  coordinate  system 
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1.   Governing  Equations  of  Motion 

The  governing  equations  of  motion  were  derived  in 
Section  II. A. 2.   They  are  restated  below. 


§|  ♦  pvv  -  o 


Dt     p 


2_<-E.)  =  o 

Dt  V 


In  polar  coordinates  with 


VA  =   ^(rA^  +  -jf] 


5_f  i  =  — (  )  +  v  —  (  )  +  —  — (  ) 

Dt1  }  3t^   '     r  3rl   '    r  301   ' 


these  equations  become 


r[p"¥  +  Pp(Vr)  +p(Vr)F]  +p(Vr)  +P(VQ)Q  +P9(VQ)   =   0     (3.24a) 


'(F[(Vr)-+  (Vr)  (Vr)_]  +  (VQ)  [(Vr)Q  -(VQ)]}   =   -r  p-        (3. 


24b) 


p{F[(VQ)-  +  (Vr)  (VQ)_]  +  (V0)  [(Ve)Q  +  (Vr)]}   =   "  PQ     (3.24c) 
-Yp{F[p"^-+  (vr)F-]  +  (vq)p"q)  +p(F[p^+  (vr)p-] 

+  (VQ)pa)   =   0  (3.24d) 
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where  p,  p,  V   and  V   are  the  density,  pressure  and  the  r- 
and  9-components  of  velocity  relative  to  the  polar  coordinate 
system  chosen  (Figure  3.10).   Bars  over  these  quantities  were 
introduced  to  indicate  their  unsteady  nature.   Quantities  t 
and  r  respresent  time  and  radial  distance  respectively  and 
bars  were  introduced  to  indicate  that  these  quantities  are 
dimensional  quantities.   Subscripts  t,  r  and  9  denote  partial 
differentiation . 

We  assume  that  the  approaching  free  stream  is  oscilla- 
tory.   Let  e  be  a  small  parameter  characteristic  of  the  departure 
of  the  free  stream  from  its  mean  (average)  constant-state 

flow.   We  also  denote  by  p,  p,  V  and  V   the  reference  steady 

r      y 

flow  quantities  (when  £  =  0)  in  the  Prandtl-Meyer  expansion 
fan  along  a  ray  9  =  constant.   We  may  express  the  flow  quan- 
tities in  the  expansion  fan  as  power  series  of  e  in  the 
following  form 

p   =  p  +  £p  +  . . .  (3.25a) 

p"  =   p  +  £p  +  . . .  (3.25b) 


V    =   V   +  £V   +  .  .  .  (3.25c) 

r      r     r 


V„   =   V   +  sV„  +  . . .  (3.25d) 


where  the  reference  steady  flow  quantities  p,  p,  V   and  V 
may  be  found  from  relations  given  in  Section  II.A.4.b. 
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We  define   non-dimensional   quantities    as    follows. 

t      =      Ujt/F,  r      =      r/T  (3.26) 

where  I   is  a  characteristic  length  taken  equal  to  one  unit 
of  length  and  U^  is  the  average  velocity  of  the  approaching 
free  stream. 

We  substitute  (3.25)  in  the  governing  equations  of 
motion  (3.24)  and  retain  only  the  zeroth  and  first-order 
terms  in  £ .   Equating  like  order  terms  in  each  equation  and 
using  (3.26)  we  obtain  the  following  two  systems  of  equations 


P(V9)Q  +  PQ(VQ)  +  p(Vr)   =   0  (3.26a) 


P(VQ) I(Vr)Q  -  (Ve)]   =   0  (3.26b) 


p  (VQ) [(Vr)  +  (VQ)Q]   =   -pQ  (3.26c) 


Y  P  PQ   =   P  PQ  (3.26d) 


r[UroPt+p(Vr)r+pr(Vr)]+p(Ve)9+p(Ve)e+p9(V9) 

+  p(Vr)+p(Vr)+pQ(VQ)   =   0  (3.27a) 


■p[uoo(vr)t+vr(vr)r]+p(ve)[(vr)9-(ve)]+P(vQ)[(vr)9 

-  (VQ)]+pV  [(Vr)Q-(Ve)]   =   -rpr         (3.27b) 


151 


rp[Uoo(V9)t+(VQ)r(Vr)]+p(VQ)[(V9)9+(Vr)] 


+  P(VQ) [(Vr)+(V9)0]+p(Ve) [(Vr)+(VQ)Q]  =  -pQ   (3.27c) 


r(p[Uoopt+pr(Vr)  ]-YPtUooPt+Pr(Vr)  ]}+(VQ)  [pp9-YPPQ] 

+  (VQ) [pp0+ppQ-Y(pp9+pp9) ]   =   0  (3.27d) 

Equations  (3.26)  constitute  the  governing  equations 
of  motion  in  the  Prandtl-Meyer  expansion  fan  for  steady  flow. 
They  are  satisfied  by  the  Prandtl-Meyer  relations  given  in 
Section  II.A.4.b  if  the  additional  relation  (VQ)  =  a  (a 
direct  consequence  of  the  fact  that  in  the  expansion  fan  the 
discontinuities  are  infinitely  weak  or  Mach  waves)  is  used, 
as  can  be  seen  by  direct  substitution.   Equations  (3.27)  con- 
stitute the  perturbation  equations  for  our  problem. 

We  assume  that  the  unsteady  parts  of  the  flow  quanti- 
ties are  of  the  following  form. 


p   =   pyMPelkt  (3.28a) 


p   =   pMReikt  (3.28b) 


VQ   =   (VQ)Ueikt  (3.28c) 


V   =   (V  )Veikt  (3.28d) 

r       r 


The  real  parts  of  the  complex  expressions  are  considered 
only. 
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where  P,  R,  U  and  Y  are  time-independent  quantities  to  be 
found  and  k  is  the  so-called  reduced  frequency  defined  by 
k  =  —  with  a)  the  circular  frequency  of  the  oscillatory 

00 

upstream  flow.   We  also  assume  that  the  amplitude  and  fre- 
quency of  oscillations  are  small  and  express  the  time-indepn- 
dent  quantities  P,  R,  U  and  Y  as  power  series  in  (ik)  of 
the  form 


P   =   p(0)  +  (ik)P(1)  +  ...  (3.29a) 


R  =   R(0)  +  (ik)R(1)  +  ...  (3.29b) 


U   =   U(0)  +  (ik)U(1)  +  ...  (3.29c) 


Y   =   Y(0)  +  (ik)Y(1)  +  ...  (3.29d) 


We  substitute  equations  (3.29)  in  equations  (3.28) 
and  the  resulting  expressions  in  the  perturbation  equations 
(3.27).   Retaining  terms  of  zeroth  and  first-order  only  and 
equating  like  order  terms  in  each  of  the  equations  we  obtain 
the  following  two  systems  of  equations. 

Zeroth  order  equations 


pM[(Vr)+(VQ)Q]R(0)  +  pM[(V9)R^0)+r(Vr)R^0)] 

+  (VQ) [p0U(O)+pU^O)]+p(Vr) [Y(0)+rY^0)]   =   0       (3.30a) 
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p(VQ)M[(Vr)e-(Ve)]R(0)+rpyMP^0)+p(V9) [ (Vr) Q- (VQ ) ] U (0) 

-p(VQ)2U(0)+p(Vr) [(VQ)Y^0)+r(Vr)Y^0)]   =   0   (3.30b) 

P(VQ)M[(Vr)+(Ve)Q]R(0)+pYMP^0)+p(VQ) [(Vr)+(VQ)Q]U(0) 

+  p(Vfl)2ui0)+rp(V J  (VQ)ul0)+p(V J  (VJY(0)  =0  (3.30c) 


(VQ)pY2MpQP(O)+ppYM[(V0)P^O)+r(Vr)P^O) ] 

-(VQ)pMpeR(0)-ppYM[(Ve)Re0)+r(Vr)R^0)] 

+  (Ve) [PPQ-YPPQ]U(0)   =   0  (3.30d) 


First  order  equations 


pM[(VQ)Q+(Vr)]R(1)+pM[(Ve)Re(1)+r(Vr)R^1)] 

+(Ve) [pU^1)+pQU(1)]+p(Vr) [Y(1)+rY^1)] 

=   -rpMUooR(0)  (3.31a) 

p(VQ)M[(Vr)Q-(VQ)  ]R(1)+rpyMP^1)+p(Ve)  [  (Vr)Q  - 
(VQ)]U(1)-p(VQ)2U(1)+p(Vr)  [(V0)Ye(1) 
+r(Vr)Y^1}]   =   -rp(Vr)UooY(0)  (3.31b) 

P  (VQ  )  M  [  (Vr )  +  (VQ  )  Q  ]  R  ( 1}  +pyMPQ(  1}  +p  (VQ  )  [  (Vr) 

+  (VQ)Q]U(1)+p(VQ)2UQ(1)+rp(Vr)  (V^U^ 

+p(Vr)  (VQ)YU)   =   -rpUoo(VQ)U(0)  (3.31c) 
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2 .  Boundary   Conditions 

The  proper  boundary  conditions  to  be  prescribed  are 
found  by  realizing  that  the  problem  is  an  initial  value  problem, 
i.e.,  the  upstream  conditions  completely  determine  the  solu- 
tion and  the  flow  matching  conditions  must  be  imposed  at  a  posi- 
tion which  differs  by  0(e)  from  that  of  OB  and  hence  must  be 
determined  as  part  of  the  solution. 

3 .  Solution 

Because  of  time  limitations  and  unsuccessful  choice  of 
test  solutions  we  have  not  been  able  to  find  a  solution. 

We  should  note  that  once  a  solution  is  found  (either 
in  closed  form  or  by  use  of  a  computer  program)  the  flow  quan- 
tities on  the  upper  surface  of  the  wedge  (Region  D)  can  be 
obtained  by  using  the  method  described  in  Section  III. A. 2.   The 
same  method  may  be  employed  to  extend  the  solution  of  the  prob- 
lem to  the  case  of  an  oscillating  wedge  in  oscillatory  flow, 
since  in  the  last  case,  the  oscillations  of  the  body  are  not 
expected  to  influence  the  flow  field  in  the  expansion  fan 
region,  except  in  the  neighborhood  of  line  OC. 
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V.   SUMMARY  AND  RECOMMENDATIONS  FOR  FUTURE  WORK 

A.   SUMMARY 

The  perturbation  method  proposed  by  Professor  Hui ,  as 
applied  to  a  two-dimensional,  oscillating  flat  plate,  with 
attached  shock  wave,  at  an  arbitrary  angle  of  attack,  in  a 
steady,  inviscid  supersonic/hypersonic  flow,  has  been  des- 
cribed in  detail.   For  periodic  oscillations  of  small  ampli- 
tude and  frequency,  the  first-order  closed  form  solutions  for 
the  flow  field  quantities  in  the  disturbed  regions,  over 
both  sides  of  the  two-dimensional  flat  plate,  have  been  given 
The  "in-pitch"  stability  derivatives  of  the  oscillating  flat 
plate  can  then  be  predicted  and  a  criterion  for  the  neutral 
damping  boundary  can  be  obtained.   The  flat  plate  results 
may  be  naturally  extended  to  include  bodies  of  slightly  more 
complicated  shapes  and  this  was  done  by  Hui  for  Caret  wings 
and  wedges  [Ref .  7] .   Utilizing  strip  theory  these  results 
may  also  be  extended  to  three-dimensional  bodies  of  a  similar 
cross  section  and  this  case  was  described  in  Section  III  for 
a  flat  three-dimensional  wing  of  arbitrary  planform  [Ref.  5]. 
The  "in-pitch"  stability  was  studied  and  the  agreement  with 
potential  flow  theory  results  (for  zero  angle  of  attack) 
was  found  to  be  good. 

The  extension  of  the  above  results  to  bodies  with  cross 
sections  composed  of  curved  segments  or  straight  segments 
forming  downstream  corners,  is  possible  (Figures  4.1  a,b) . 
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In  Section  IV  the  perturbation  method  proposed  by  Hui 
was  extended  to  include  the  effects  of  periodically  oscillating 
upstream  flow  over  both  sides  of  a  stationary  flat  plate. 
The  method  may  be  readily  extended  to  include  oscillations 
of  the  plate.   The  assumed  form  of  upstream  flow  oscillations 
includes  as  special  cases  the  form  of  induced  oscillations 
in  the  flow  field  by  the  oscillating  flat  plate  considered 
before  in  steady  supersonic/hypersonic  flow.   The  solution 
over  the  compression  side  is  given  in  closed  form,  while 
that  over  the  expansion  side  is  given  as  a  series  of  iterative 
calculations . 

B.   RECOMMENDATIONS  FOR  FUTURE  WORK 

The  upstream  oscillatory  flow  solutions  may  be  combined 
with  the  upstream  steady  flow  solutions  to  study  three 
general  problems. 

First,  the  induced  flow  field  by  two-dimensional  bodies, 
of  any  cross  section,  at  arbitrary  angles  of  attack  and 
attached  shock  waves,  in  steady  flow.   For  this  problem  the 
flow  field  can  be  obtained  in  steps,  from  leading  to  trailing 
edge,  over  the  straight  line  cross  section  segments  (Figure 
4.1a).   Curved  segments  could  also  be  approximated  by  a 
series  of  straight  line  segments  (Figure  4.1b). 

Second,  the  induced  flow  field  by  a  formation  of  bodies 
with  the  exception  of 

i)   Regions  behind  crossing  shocks, 
ii)   Expansion  fan  overlapping  regions, 
iii)   Mixed  regions  of  above. 

157 


^ 


shock 


(a) 


REGIONS     B 


(  b) 


Figure    4.1.      Illustration   of   flow  solutions 
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Third,  the  induced  flow  field  behind  three-dimensional 
bodies  of  any  cross  section,  utilizing  the  strip  theory- 
approximation  . 

The  dynamic  stability,  in  pitch,  of  the  bodies  involved 
in  each  of  the  above  general  problems,  might  be  studied. 

A  closed  form  solution  for  the  expansion  fan  problem 
formulated  in  Section  IV. C  might  be  sought.   Similarly 
a  closed  form  solution  might  be  sought  for  the  more  general 
expansion  side  problem  considered  in  Section  IV. B.   This 
solution  might  possibly  be  found  by  letting  the  step  angle 
9,  (representing  the  step  change  in  flow  direction)  tend 
to  zero  in  the  basic  formulas  that  are  iteratively  used. 
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APPENDIX  A 
PROCEDURE  FOR  SOLVING  SHOCK  BOUNDARY  CONDITIONS 


Let  Z1  ■  I   gradG 

X   =  J[^-   +  V   •  gradG] 

1     at    X 

X9  =  ![—  +  V,  •  gradG] 

*  3t     Z 

where  I    is  the  characteristic  length  chosen,  G  =  0  is  the 
equation  of  the  shock  and  the  subscripts  1  and  2  denote 

quantities  before  and  after  the  shock  respectively. 

2  Y 
Equations  (2.3a-2.3c)  are  written  with  y  =  — — j- 


P1X1   =   P2X2  (A,1) 

2       2         2       2 
p  X   +  p  z    =   P2X2  +  P2Z  (A. 2) 

2      ^1    2        2      P?   2 
X*  +  y  -±.  zz   =   x;  +  y  -=•  z^  (A. 3) 

1      pl  2      p2 

Consider  a  harmonic  oscillation  with  frequency  co  and  maximum 
angular  deviation  of  the  unsteady  flow  from  the  steady  reference 

flOW  £ . 

Let  e  =  ee    where 

u 
t  (non-dimensional  time)   =  —  t 


k  (reduced  frequency)      =   — 
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t:   dimensional  time 
u  :   reference  velocity  (in  steady  flow) 

For  the  general  case  assume  that  the  approaching  flow  is  also 
unsteady  and  let 


2 
Xl   ~   Xll  +  £X12      X2   "   X21  +  £X22    Z    "   Zll+eZ12 


Pl   =   Pll  +  £p12      P2   =   P21  +  £P22 


pl   =   pll  +  £p12      p2   =   p21  +  £p22 


We  substitute  these  values  in  (A.1-A.3)  and  retain  only  zeroth 
and  first  order  terms  in  e .   Equating  zeroth  and  first  order 
terms  we  get  the  following  two  systems  of  equations  respectively 


P21X21   =   P11X11  (A. 1.1) 


P21  [P2i+Pii(1~y)  J   =   Pn  [pii+p21(1""y)  1  (A. 2.1) 

2  2   2 

P21tpllXll+ypllZll1   =   (U-1)P11X11  (A. 3.1) 


p21X22   =   P11X12  +  P12X11  "  P22X21  (A'1-2) 


p22  (P21-kpxl)   =   p12  (P11+kp21)  -P22  (Puk+P21)  +Pi2  (pll-kp21} 

(A. 2. 2) 
P22(pllXll+ypllZll)   =   2pllXllX12[(y-i:'pirp21] 


+  P12Xi1[2(y-l)p11-P21]-pp21[p11Z12+p12Z11]  (A. 3. 2) 


with  k  =  p-1  =  (y+1)/(Y-D. 
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Equations  (A. 1 .1-A. 3 .1)  give  the  unknown  flow  quantities 
behind  the  discontinuity,  for  the  steady  reference  flow  (in 
terms  of  known  flow  quantities  before  the  discontinuity)  and 
are  also  used  to  simplify  equations  (A. 1 . 2-A. 3 . 2)  even  further 
These  quantities  will  then  be  used  in  simplified  equations 
(A.l . 2-A. 3 . 2)  to  give  the  unknown  flow  quantities  behind  the 
discontinuity  for  the  unsteady  flow. 

To  complete  the  system  of  boundary  conditions  across  the 
shock  we  now  consider  equation  (2.3d)  which  expresses  the 
conservation  of  tangential  momentum.   This  equation  is 
equivalent  to 


Vl  *  T   =   V2  "  T  (A. 4) 


where  T  is  a  tangent  to  the  surface  vector,  such  that 
VG-T  =  0. 

For  a  surface  G  =  x  tan  4>  -y  +  eQ(x)  =  0  we  will  choose 
this  vector  to  be  T  =  {l,tan  <j>  +eQ'}. 
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APPENDIX  B 
OSCILLATING  WEDGE  IN  STEADY  UNIFORM  FLOW 

For  this  case  which  was  considered  in  Section  III.A.l 
(Figure  2.1a) , 


V,   =   {u,,v,}   =  U^icos  0,-sin  6} 


V0   =   {u0/v0}   =  u  {1  +  eU,  eV} 

2         2   2         O 


G   =  x  tan  $  -  y  +  eQ  (x) 


Then 


VG   =  -{tancf)  +£0.-1),  T  =   {1,  tan  $  +eQ  > 

I 

2  2 

Z   =   Z11  +  eZ12   =   1  +  tan  <j>  +  e[2  tan$Qx] 

X,   =  X„„'  +  eX,„   =  U  {sin  I  +  e[cos  9Q   +  ~(ik)Q]> 


'11     12      °°  cos  <j>  x   U 


OJ 


X2   =   X21  +  £X22   =  u0{tan  *  +  £  [U  tan<^  "V  +Qx  +  (ik)Q^ 


We  also  have 

Pl      =      Pll    +    £p12      =      P~   +   £*°        P2      =      p21+£p22=Po+£poYMoP 
pl      "      pll   +    £p12      =      p~   +   £'0         p2      =      p21+£p22=po+£poMoR 


163 


These  values  will  be  substituted  in  equations  (A.1-A.3) . 

Equations  (A.l .1-A. 3 .1)  are  considered  first.   From  Equation 

2 
(A.  3.1)  we  get  after  multiplying  it  by  cos a  j- 

P    U 

^^^  oo     oo 
po  2  2        Y     Poo  Y+l 

-£[sin^3   +  -~  y]      =     £J  sin^B  (A  1.1) 

00  Y         p    U  ' 

00       00 

Y     P-  1  PQ     • 

Setting  «-  =  — y  and  solving  for  —  we  get 

p  U     M  P« 

r»  CO        oo 


P  (y+DM2sin23 

O  00 


Poo      (Y-l)sin23M2  +2 


(A.3.1.A1) 


From  equation  (A. 2.1)  we  get  after  dividing  it  by  p   and 
solving  for  — 

t"^oo 


?2    -     n   -  ^£  X±Ii/r^°  -  I±Ii 

PM  P   Y-l  7  P     Y-l 

roo  ^  00   '         ^  oo      ' 

Substituting  for  p  /p   its  value  from  (A.3.1.A1),  we  get 

2   2 
p      2Ysin  3M  -  (y-1) 

(A.2.1.A1) 


Pw  Y+l 

From  equation  (A.  1.1)  we  get  after  multiplying  it  by  cos  <J> 


p  U  sin  6   =   P  u  sin  d>  (A  1.2) 

00   CO  OO 


From  equation  (A. 4)  we  get  the  following  two  equations  equating 
the  zeroth  and  first  order  terms  in  e  and  using  the  relation 
cos  9  cos  tj)  -  sin  9  sin  <j>  =  cos  3 


U  cos  6   =   u   cos  <b  (A. 4.1) 

00  o 
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-U   sin  6Q        =      u    (U+tan<J>V)  (A. 4. 2) 

From  equations    (Al.1.0)    and    (A4.1)    we   have 

TT     '     HH  (A.1.1.A1) 

p  tan   <p 

00 

Equations  (A.1.1.A1),  (A.2.1.A1)  and  (A.3.1.A1)  consti- 
tute the  supersonic  uniform  wedge  flow  solution  presented  in 
Section  II. A. 4. a. 

Equations  (A.l . 2-A. 3 . 2)  are  considered  next.   The  time  - 
independent  unknown  functions  R,  P,  U  and  V  can  be  found  in 
terms  of  Q1  and  (ik)Q  in  the  following  order.   From  simplified 
equation  (A. 3. 2)  the  function  R;  from  the  known  function  R 
and  equation  (A. 2. 2)  the  function  P;  finally  from  equations 
(A.l. 2)  and  (A. 4. 2)  the  functions  U  and  V.   The  derivations 
are  lengthy  and  tedious  and  will  not  be  included  here.   The  ex- 
pressions for  the  complete  set  of  coefficients  in  equations 
(2.11)  are  given  in  [Ref.  7]  and  are  quoted  below. 

(l-W)A     =      (1 — -)cos%[l+—  (M  -l)tan^<j>-YW(—  -  1)  ]  (A1.3a) 

P0  w  oo 

p  0  p  p 

(l-w)B      =       (1  -~)cos    <|>[1  +~  W-yW(— -1)]  (A1.3b) 

^O  °o  co 

Pco  4  Y-l         Pn 

(l-W)C      =      2kH(l )cos%[l  -2~W(-2.  -1)  ]  (Al.4a) 

po  p°o 

D      =      C/cos2<})  (A1.4b) 

E      =       (A-B)    cot   d>  (A1.5a) 
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F   =   -B  tan  <j>  (A1.5b) 


p      4  po 

(l-W)G   =  kH(l  --2.)cos%[y+1-(Y-1)-°-]  (A1.6a) 

P0  P°o 

J  =   G/cos2<j>  (A1.6b) 

M  

k  =  —      ,   H  =   VM  -1  tan  <J>,   W  =  M  sin  1 

o        r  o     Y 


n/m2  -1 


o 

To  illustrate  the  derivation  procedure  suggested,  we 
find,  below,  the  expression  for  the  function  R  and  compare  the 
coefficients   G  and  J  with  those  given  by  relations  (A1.6a) 

and  (A1.6b)  above. 

2     2 
We  multiply  equation  (A. 3. 2)  by  cos  ^/P^U^  and  use  relations 

(A.1.4.A1),  (A1.2),  (A. 4.1),  (A.1.1.A1)  and  cos  9  =  cos  (g-<f>) 

to  simplify  it.   We  finally  have 

o         o    ' 
Thus  the  coefficients  G  and  J  found  are 

~G     =  l   «o»V  (!=.  -  IZl)  (Al.7a) 

M  sm  0  p     Y+l 
o     T   o    ' 

J  =  G/cos24)  (Al.7b) 

The  second  coefficient  J  is  expressed  by  the  identical  rela- 
tions (A1.6b)  and  (A1.7b).   It  can  also  be  seen  that  the 
expressions  (A1.6a)  and  (Al.7a),  for  the  first  coefficient 
G,  are  identical  by  direct  substitution  and  use  of  the 
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following  additional  relations  given  in  Section  II. A. 4. a, 


poo     2  +  (y-Dx1  2  +  (Y-l)x' 


where 


=  M  sin  8/ 


2    2 

x   =  M  sin  <p 
o    r 


x   = 


PQ       (Y+Dx'     '      ~      2yx'  -  (Y-D 
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APPENDIX  C 

OSCILLATING  FLAT  PLATE  IN  STEADY  UNIFORM  FLOW — 

EXPANSION  SIDE 

For  this  case,  which  was  considered  in  Section  III. A. 2 
(Figures  2.2a  and  2.2b)  we  assume  that  there  is  no  change 
in  the  steady  flow  quantities  across  the  front.   Thus 


Vl   =   ui{1'0}  '     v2   =  ul{1  +eU'eV} 


G   =   x  tan  <|>  -  y  +  eQ  (x) 


Then 


VG     =     -{tan   <J>   +  eQ',-1},  T     =      {l,tan   <J>   +   eQ '  } 


2  2 

Z   =  Z11   +   eZ12   =   1  +  tan  <J>  +  e[2  tan^Q'] 


Xl   =   Xll  +  £X12   =  ui^tan  *  +  etQ'  +  (ik)Q]> 


X2   =   X21  +  £X22   =  ui^tan  ^  +  e[U  tan  <J>  -  V  +  Q '  +  (ik)Q]> 


We  also  have 


Pl   =   Pll  +  £p12   =   Pl  +  £'° 


P2   =   P21  +  £P22   =   Pl  +  £pl^MlP 


Pl   =   Pll  +  £p12   =   Pl  +  £ "  ° 


168 


P2   =   P21  +  £p22   =   pl  +  £e]JMlR 


Equation  (A. 3. 2)  is  simplified  using  equation  (A. 3.1)  and 
becomes  with  p.-  =  0, 

»22    -    pai«i-?ir^>«^-l5f» 

Subsituting  values  given  above  we  get 


R      =      G'Q'    +   J' (ik)Q  (A2.1) 

with 

G'    =    rFifspk*  '   J'   -  =vcos2*  (A2.1.D 

Equation  (A. 2. 2)  is  simplified  using  equation  (A. 2.1) 
and  beomces  with  p.  2  =  p.  2  =  0 

pn  (PU^P12) 
P22      p22  pxl  (p2]_k  -Pll) 

We  substitute  values  from  above  and  have 

P   =   R  (A2.2) 

From  equations  (A. 1.2)  and  (A. 4)  we  get  the  following  two 
relations  respectively, 


U   tan   <J>    -  V      =      -M     tan  <j>    R 


U  +  V  tan    c|>   =      0 
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Solving   these   relations  we   get 

V      =      A'Q'    +   B1 (ik)Q  (A2.3) 


with 


4 
A,      =      4cos    (j)    ^  B,      =      A,/cog2(J)  (A2.3.1) 

T  "**■■ 


and   also 


U      =      E'Q1    +  F1 (ik)Q  (A2.4) 


with 


3 

t?i                  4cos    djsin  d>  „ ,  _,  ,        2, 

E*       = p~£f} *-    /         F'      =      E'/cos    <J> 


(A2.4.1) 


Equations  (A2.1-A2.4)  give  the  boundary  conditions 
(2 . 16f-2 ,16i)  and  (2 . 17f-2 .17i)  when  the  power  series  expan- 
sions of  the  time -independent  quantities,  given  by  relations 
(2.8),  are  used.   The  apparent  differences  between  the  ex- 
pressions given  here  for  the  boundary  conditions  and  those 
given  in  [Ref .  8] ,  formulas  (21)  and  (22)  ,  are  due  to  the 
different  representations  of  the  equation  of  the  front. 
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APPENDIX  D 
OSCILLATORY  UPSTREAM  FLOW  PAST  A  STATIONARY  WEDGE 

For  this  case  which  was  considered  in  Section  IV. A. 3 
(Figure  3.4b)  we  obtain  from  relations  (3.6)  and  (3.9)  with 
£2  =  0, 


vi    =    v{cos  9d  +  £iV-sin  9d  +  £lV 

B 
Pl      =      P11+£1P12      =      Pon+£lPouPB  (A3'la) 

pl      =      pll+£p12      =      poD+£lpo*B  (A3*lb) 

B  B 

v2     =     uo{l+£lU1,e1V1} 


p2     =     P21+e1P22     =     P0+^P0YMoPl  (A3*lc) 

p2      =      P21+£1P22      =      P0+^iP0MoRl  (A3. Id) 

where  U    ,    V    ,    P„   and   Pw   have   the   following   general  form  given 
by    (3.6a-3.6d) 


M_      =      C       +  (ik,  )  (C      x  +C      y  +C      )  (A3.0) 

B  mQ  1        m1  m2  m3 


We 


also  have  with    %2   =   I    and  G   =  x  tan  <f>  -y  +  e,Q,  (x)    =   0 

VG      =       ~{tan  tj>  +  £^{,-1},         T      =      { 1 ,  tan  <J>  +   E^Qp 
I 

Z2      =      zn+£izi2      =      1  +  tan2(J)  +e1*2  tanc})Q|  (A3.1e) 
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X  =     u      sin   3/cos    (J>  (A3. If) 

°B 

X12      =      uo    [UB  tan^-vB  +cos    9dQl+(ikl)uoQl/uo    ]  (A3-1(?) 

B  B 

X21      =      Uo  tancJ)  (A3,  lb) 

X22      =      UotUl  tan<^  "  Vl  +Q1  +  (iki)Qil  (A3.1i) 

Equation    (A. 3. 2)    is    simplified  using   equation    (A. 3.1) 
and  becomes 

n       r      X12      p12      Z121M         1      p21n    ,  p12p21      yp21P12Zll 

p    0    =    P01  [2   ^ —  +  - - — J  [1  -  — T J   + — = - 

22  21        Xu      pn      Z1X  u-1    pn  pu  p^^ 

Using   relations    (A3.1)    we   obtain 


Rl      '      KR1UB+KR,VB+KR,RB+VB+VQi+KRfi(ik>Ql  (A3-2a> 

12  3  4  5  6 

with 


K_,        =     H,  sin  <t>/sin   3,        K_.        =      -H,   cos  4>/sin   3, 
Rl  l  R2  L 


K  =      1/M     +H,/2,  K  = ^ 2 

K3  o        l  K4  M    [sin^3(Y-l)Mz     +2] 


°  °B 


Kn        =     H,  cos2<J>/tan   3    ,      K_        =     H,/tan   3 
R5  1  R6  l 


y-1    po 

Hi  -  2(1  -Y^r)/Mo 

°B 
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Equation  (A. 2. 2)  is  simplified  using  equation  (A. 2.1) 
and  becomes 

p21  [Pn+(^"1)P21]       P21  [Pn+(^l"1)P2l] 
P22   "   p12  p^  LPll-(u-l)p21]  "  P22  p~^  (PU-(H)P21] 


+  P 


P_21 
12  pn 


We   substitute   values    from    (A3.1)    and   get 


P.       =      K      U+K      V„  +  K      R„  +K      P„  +K      Q'+K       (ik.  )Qn        (A3. 2b) 
1  p]_    B         p2    B         p3    B         p4    B         p5vl         p6         1    vl 

where : 


Kp         =      H2KR_     ,         j    «    1,2,5,6 

i  +  (u-Dp0/p0b 
KP3     =     H2KR3  +     Mo[l-(y-l)P0/POBJ 

Kp4      =      H2KR4    +    l/(YM0)f 

p0[(y-i)+p0/p0Bl 

H2      YPn  [1  -  (y-l)p n/Pn  J 

°B  B 

From  equations  (A.  1.2)  and  (A. 4)  we  get  a  system  of  two  equa- 
tions and  solving  for  U,  and  V.  we  obtain 


V,       =      K      U+K      V+K      L+LP.+K     Q.'+K     Q.liL)  (A3. 2c) 

1  v,    B      v2   B      v^    a     " .    B      v5    1      Vg    1         1 

where 
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V       =      H3KR.        ^    "    3'4 


K  =      H,K_      +   H„,               K           =      H,K_      +    H./tan    9, 

vl  3    Rl           4                   v2                 3    R2           4                d 

Kv  =      H3KR      +  H4    cos(3+<J>)/tan   cj> 

5  5 

K  =      H.K        -    cos2<Hpo   /p     -1] 

v6  3    R6                          °B      ° 

H_  =     M     tan  (J) ,             H.      =      sin   9,  sin  <$>  cos  <£ /sin    3  cos     3 

3  o                             4                         d 


Ux      =      Ku  UB+KU   VB+KU   Rg+Ku   PB+KU   Q^+Ku   Q1(ik)  (A3 . 2d) 

with 


Ku        =      tan  cj>  Kv         (j    =    3,4,6) 


K  =     K       tan<£  +u      /u    ,  K  (-K       tan  0  +  u     /u   )  tan   0 

ul  vl  °B      °  U2  v2  °B      ° 


K  =      -(K       tan(})+u       sin9,/u    ) 

u5  v5  oB  d'    o 


The  general  form  of  the  shock  boundary  conditions  expressed 
by  equations  (A3. 2)  is 


Yl   =   KY1VKY2VKy3RB+KY4PB+KY5Qi+KY6Ql<ikl)         (A3'3) 


Using  relation  (3.10)  the  zeroth  and  first  order  boundary 
conditions  across  the  shock  become 


Y<°>   =  K^u'O^v'^K^r'^p'0^;'0'      (A3.4a) 
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'i11  -  VB1)+VB1S^>+v*1)+via,+Vi0) 


where 


(A3. 4b) 


Introducing    (A3.0)    we   obtain 


mJ0)       =      C  and        MTf1)       =      C      x+C      y+C 

o  1  2  J 


Relation  (A3. 4a)  then  gives 


y|0)   =   K   Q[(0)  +  K     At  y  =  x  tan  0  (A3. 5a) 

x 


Kv         =      K      C         +    Kxr    C         +    K_,    C         +    K      C 
Yx  *1    Uo  Y2    Vo         ^3    ro  Y4    ?o 


Similarly,    relation    (A3. 4b)    gives 

Y^1]       =      Ky      x+I^      y  +Ky  +  Ky   Q^  (1)   +  K^   Q  (A3. 5b) 

xx  yy  x  o 

At     y   =   x   tan   <j> 

where : 


Kv  =      KTr   C        +   Kv   C        +   K__    C        +   K_,   C 

Yxx  Yl   Ul  Y2   Vl        ^3    rl        %    Pi 


VL,  =      IL   C         +    Kv    C         +    K^   C         +    Kv    C 

Y  Y,    u_  Y~    v0  Y0    r0  Y,    p0 

yy  12  22  32  4^2 
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Kv  ■  Ky  Cu  +  KY  Cu   +  KY  Cr   +  KY  Cp 
Y      Yl  u3     Y2  u3     Y3  r3     *4  p3 


Ky    =   K    ,       K     -   K 
x        5  ob 
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APPENDIX  E 
IN-PITCH  DERIVATIVES  BASED  ON  TN ' S  2699  AND  319  6 

To  find  the  stiffness  and  damping-in-pitch  derivatives 
given  in  [Ref.  5  ]  we  proceed  as  follows. 
We  define  the  following  quantities: 

S  =  planform  area 

b  =  half  wing  span 

C  =  root  chord 

C,  =  tip  chord 

L_E  =  leading  edge  angle  of  sweep 

M  =   free  stream  Mach  number 

2    0  5 

3   =   cotangent  of  Mach  angel  =  (M  -1) 

2 
A  =   aspect  ratio  =  (2b)  /S 


A 


=   taper  ratio  =  C  /C 


_  2 

C   =  mean  aerodynamic  chord  =  2C  (A  +A+1)/3(A+1) 

m  =   slope  of  LE  =  cot  L 

Li£ 

A*   =  parameter  =  A3 
m'   =  parameter  =  m3 

We   also   state   the   following   stability  derivatives   defined 
explicitly   in    [Ref.    11,12] . 

Resulting   from  steady   state   motion 


C 
L        =      lift  curve-slope   derivative 

a 

CT  ,  C    =   lift  and  pitching  moment  derivatives  (due  to 
q   mq  steady  pitching) 
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Resulting  from  time-dependent  motion, 


C_    =   lift  derivative 

a 


C  .   =  pitching  moment  derivative  (due  to  constant 
a     vertical  acceleration) 


To  account  for  differences  in  definitions  of  quantities 
we  form  the  following  factors 


Fx   =   C/Cr,  F2   =   F^/2 


The   stiffness    and  damping-in-pi tch  derivatives    as    defined 
in   this   thesis    are    then   given  by 


C  =      F,C         ,  C  =      F0(C       +C      )  (Bl) 

mQ  1  ma  m§  2     mq        m^ 

To  find  C   ,  C   ,  C    from  NACA  TN's  we  proceed  as  follows 
m    m    m  • 
a    q    a 


i)   Form  3 


,  A  ,  A1 ,  m1 ,  cot  m* ; 


ii)   Read  from  NACA  TN  2699  (Figures  11-15)  3C   and  find 

m 
a 

C    =  (j3C   )/S 
m       m 
a       a 

iii)   Read  from  NACA  TN  2699  (Figures  26-30  3CT   and  use 

Li 

a 
its  value  to  find  pivot  position  as  follows 


X/C   =   -F. (SC   )/(3C.  ) 
r       1   m      L 
a      a 

iv)   Read  from  NACA  TN  2699  (Figures  16-20  and  Figures  21-25) 

3C*  and  3C*  which  refer  to  pivot  position  found  above. 
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v)   Use  following  formulas  to  find  3CT   and  £C 

Lq       "q 

°L   =   (6CL   "  26Cm  )/S 
q  q        a 


C    =   [BC*   +  (BC   )(BCT  )/(BCT  )]/B 
ni         m       m     L      L 

q  q        a      q       a 

vi)    Read  from  NACA  TN  3196  (Figures  11-15)  quantities 

($C   )   and  ($C   )^  and  use  them  to  find  C   as 
m*  i        m •  2  m* 

a  ±  a  a 

follows, 


C    =   [M2(SC   )-/B2  +  (M2/S2  +1)  (6C.)9]/B 
m*  m  •  l  m  •  z 

a  a 

a 

vii)   Use  relations  Bl  to  find  comparison  values  for  C 

m( 

and  C 

me 
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